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test  methods  which  are  in  use  are  described  in  MIL-STD-883B.  The  present  study 
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microelectronic  package  during  a thermal  shock  test.  In  this  model,  an 
approximate  model  is  adopted  where  the  lead,  the  glass  seal  and  a part  of 
the  base  are  represented  by  3 concentric  cylinders. 

In  order  to  carry  out  the  thermal  stress  analysis  a temperature 
analysis  is  performed  to  determine  the  unfiorm  temperature  of  the  lead, 
the  uniform  temperature  of  the  base  and  the  radial  temperature  distribution 
in  the  seal.  This  temperature  model  takes  into  account  the  entire  amount 
of  the  base  material  and  lead  material  in  the  package  and  therefore  is  some- 
what different  than  the  thermal  stress  model. 

The  temperature  model  can  be  used  to  determine  the  approximate 
temperature  distribution  in  the  three  cylinder  model  if  the  effective  heat 
transfer  coefficient  that  determines  the  rate  of  heat  transfer  between  the 
lead  and  the  fluid  and  between  the  base  and  the  fluid  is  known.  The 
effective  heat  transfer  coefficient  is  evaluated  using  experimentally 
obtained  time-temperature  histories  for  the  lead  and  the  base. 

The  analysis  is  applied  to  a commercially  used  microelectronic 

package.  The  transient  thermal  stresses  developed  in  the  seal  during  the 

thermal  shock  test  level  A of  MIL-STD-883B  are  computed  and  presented  in 

the  form  of  graphs.  The  analysis  is  then  applied  to  packages  with  different 

geometries.  The  results  thus  obtained  are  non-dimensionalized  and  presented 

in  the  form  of  graphs.  These  graphs  can  be  used  to  determine  the  thermal 

stresses  in  any  package  for  which  the  ratio  of  the  material  properties 

(Young's  modulus,  coefficient  of  thermal  expansion,  Poisson's  ratio,  heat 

transfer  property  ratios)  of  the  sealing  material  to  those  of  the  metal 

are  equal  to  the  corresponding  material  property  ratios 

((hL(A./N)A  (L-L  )2  = 96.459,  p c /p  c = 1.681,  E /E  = 2.439, 
d se  o m m se  se  m se 
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1.269,  v /v  = 1.364)  of  the  experimental  package  used  in  this 


investigation.  The  present  study  is  of  general  use  to  designers  and  screeners 
of  microelectronic  packages  since  the  same  combination  of  materials  as  the 
ones  of  the  experimental  package  is  used  in  most  packages. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOEfWlAn  Dm*  Emma*; 


I 


iii. 


ACKNOWLEDGEMENTS 


The  present  report  is  based  on  the  thesis  submitted  by 
the  first  author  to  the  Graduate  School  of  Syracuse  University 
in  partial  fulfillment  of  the  requirements  for  the  degree  of 
Master  of  Science  in  Mechanical  Engineering. 

Financial  support  for  this  research  was  provided  by  the 
Rome  Air  Development  Center  under  Grant  Number  F30602-78-C-0083 , 
Task  Number  PRN-8-5165. 

The  authors  wish  to  express  particular  appreciation  to 
Robert  W.  Thomas  of  RADC  for  identifying  the  problem  that 
initiated  the  study  and  also  to  Peter  Manno  of  RADC  for  his 
helpful  discussions. 

The  authors  would  also  like  to  thank  the  Isotronics  and 


Tekform  microelectronic  package  manufacturing  companies  for 
donating  the  packages  used  during  experiments,  the  Bendix 
Electrical  Components  Division  for  contributing  the  pressure 
bomb,  and  J.  C.  Gioia  from  General  Electric  Company  for  his 
valuable  suggestions  concerning  the  sealing  of  the  experimental 
packages. 

Special  thanks  to  Robert  P.  McEvoy  for  helping  during  the 
writing  of  the  manuscript  and  to  Cynthia  Porter  and  Kathy  Curtice 


ACKNOWLEDGEMENTS i 

ABSTRACT iii 

SYMBOLS  AND  UNITS  v 

CHAPTER  I : INTRODUCTION 1 

CHAPTER  II  : THERMAL  STRESS  ANALYSIS  10 

CHAPTER  III  : TEMPERATURE  ANALYSIS 16 

CHAPTER  IV  : EXPERIMENTS 23 

4.1  Procedure 


4.2  Evaluation  of  the  Effective  Heat 
Transfer  Coefficient 

4.3  Experimental  Results 

CHAPTER  V : RESULTS 48 

CHAPTER  VI  : PRACTICAL  APPLICATION  OF  THE  RESULTS.  . 49 

CHAPTER  VII  : CONCLUSIONS  AND  RECOMMENDATIONS  ....  85 


APPENDIX  A : FINITE  ELEMENT  ANALYSIS  FOR  EVALUATING 

THE  TRANSIENT  TEMPERATURE  DISTRIBUTION 
IN  THE  GLASS  AND  RELATED  APL  PROGRAMS 


A . 1 Evaluation  of  the  Transient 
Temperature  Distribution  in 
the  Glass  Using  the  Finite 

Element  Analysis  91 

A. 2 Programs  "TEMPDIS"  and  "FEM" 
to  Determine  the  [K] , [C]  and 

[F]  Matrices 95 

A. 3 Program  "GTS"  to  Compute  the 
Transient  Temperature  Distri- 
bution in  the  Glass 96 


APPENDIX  B : APL  PROGRAMS  COMPUTING  THE  TRANSIENT 

THERMAL  STRESSES 

B.l  Main  Program  "STR"  to  Compute 
the  Transient  Thermal  Stresses 


at  a Given  Radius 98 

B.2  Programs  "STRESS"  and  "Z"  to 
Compute  the  Stresses  at  each 
Increment  of  Time 100 


APPENDIX  C : APL  FUNCTION  EXPLAINING  HOW  THE  THERMAL 

STRESSES  ARE  COMPUTED  USING  THE  APL 
PROGRAMS  PRESENTED  IN  APPENDICES  A 
AND  B 

C.l  Function  that  Explains  the  Use  of 
STR  and  Illustrates  the  Use  of 
the  Program  with  an  Example.  . . .102 

REFERENCES 105 


1 
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SYMBOL 


DESCRIPTION 


Thermal  diffusivity  of  glass  (k  /p  c ) in  /sec 

g g g 

Outside  surface  area  of  the  lead  emerging  from 
the  base  in2 

2 

Surface  area  of  glass-base  interface  in 

2 

Surface  area  of  lead-glass  interface  in 

2 

Outside  surface  area  of  the  base  and  the  lid  in 

Specific  heat  of  chromel  thermocouple  wire 
Btu/lb- °F 

Specific  heat  of  constantan  thermocouple  wire 
Btu/lb- °F 

Specific  heat  of  glass  Btu/lb-°F 
Specific  heat  of  Kovar  Btu/lb- °F 
Specific  heat  of  any  metal  Btu/lb-°F 

Specific  heat  of  solder  Btu/lb- °F 

Specific  heat  of  any  sealing  material  Btu/lb°F 

Diameter  of  lead  in 

Diameter  of  chromel  thermocouple  wire  in 

Diameter  of  constantan  thermocouple  in 

Diameter  of  solder  sphere  in 

Young's  modulus  psi 

Young's  modulus  for  glass  psi 

Young's  modulus  for  Kovar  psi 

Young's  modulus  for  any  metal  psi 

Young's  modulus  for  any  sealing  material  psi 

Non-dimensional  Young's  modulus  ratio  (E./E  ) 

k g 

n 2 

Least  square  error  defined  by  .1,  (y.-y  .) 

i=l  l ei 

2 

Effective  heat  transfer  coefficient  Btu/hr-ft  -°F 
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Non-dimensional  effective  heat  transfer  coefficient, 
K*o 


Thermal  conductivity  of  glass  Btu/sec-in-°F 


Thermal  conductivity  of  any  sealing  material 
Btu/sec-in-°F 

Total  length  of  the  lead  in 

Length  of  the  lead  emerging  from  the  base  in 
Length  of  the  thermocouple  wires  in 
Total  number  of  leads 
Radial  coordinate  in 

Outside  radius  of  the  third  outside  cylinder  in  the 
3-cylinder  model  in 

Outside  radius  of  the  glass  annulus  in 
Radius  of  the  circular  lead  in 

hA 

Non-dimensional  radial  parameter,  - — 2—  r 

kgA2 

Temperature  °F 

Initial  uniform  temperature  of  the  package  °F 

Temperature  of  the  fluid  °F 

Displacement  in  the  radial  direction  in 

3 

Volume  of  the  base  and  lid  in 

3 

Volume  of  the  chromel  thermocouple  wire  in 

3 

Volume  of  the  cons tan tan  thermocouple  wire  in 

3 

Volume  of  the  lead  in 

3 

Volume  of  the  solder  in 

Any  temperature  on  the  theoretical  time- 
temperature  curve  °F 

Any  temperature  on  the  experimental  time- 
temperature  curve  °F 
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Coefficient  of  thermal  expansion  in/in-°F 
Coefficient  of  thermal  expansion  of  glass  in/in-°F 

Coefficient  of  thermal  expansion  of  Kovar  in/in-°F 


Coefficient  of  thermal  expansion  of  any  metal 
in/in- °F 

Coefficient  of  thermal  expansion  of  sealing 
material  in/in-°F 

Non-dimensional  coefficient  of  thermal  expansion 
ratio,  am/ase 

Strain  in  the  radial  direction  in/in 

Strain  in  the  tangential  direction  in/in 

Strain  in  the  axial  direction  in/in 

Difference  between  final  and  initial  temperature 
during  any  temperature  excursion  °F 

Difference  between  final  and  initial  temperature 
during  a thermal  shock  test  °F 

Difference  between  room  temperature  and  the  lower 
bath  temperature  of  a thermal  shock  test  °F 

0 

Non-dimensional  temperature  difference,  — 

o 

Poisson's  ratio 

Poisson's  ratio  of  glass 

Poisson's  ratio  of  Kovar 

Poisson's  ratio  of  any  metal 

Poisson's  ratio  of  sealing  material 

Non-dimensional  ratio  of  Poisson's  ratios,  v /v 

m se 

3 

Density  of  glass  lb/in 


viii 


Density  of  the  chromel  thermocouple  wire  lb/in3 
Density  of  the  constantan  thermocouple  wire  lb/in3 
Density  of  Kovar  lb/in3 

Density  of  any  metal  lb/in3 

. 3 

Density  of  solder  lb/in 

Density  of  any  sealing  material  lb/in 

Thermal  stress  in  the  radial  direction  psi 
Thermal  stress  in  the  tangential  direction  psi 


o Thermal  stress  in  the  axial  direction  psi 

z 

t Time  sec 

t Reference  time  arbitrarily  chosen  to  be  the  time 

r necessary  to  bring  to  90%  of  its  equilibrium 

temperature  1 lb  of  Kovar  when  it  is  heated 
from  32 °F  to  212°F  with  h = 210  Btu/f t2-hr-°F 

hAo 

t*  Non-dimensional  time,  t j-  t 

ckpkV£ 


EVALUATION 
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\ 

The  objective  of  this  effort  was  to  develop  a mathematical  model 
for  predicting  the  transient  thermal  stresses  that  occur  in  a typical 
annular  lead-through  glass  seal  of  a microelectronic  package  during  a 
thermal  shock  test. 

The  study  result?  show  that  the  magnitude  of  the  thermal  stresses  in 
the  seal  area  depends  upon  the  initial  residual  stresses  that  result  from 
cooling  the  package  to  room  conditions  after  the  high  temperature  sealing 
process  is  completed.  A computer  program  was  developed  to  analyze  the 
thermal  stresses  in  any  circular  lead  package  while  undergoing  a thermal 
shock  test.  A method  has  been  developed  for  predicting  the  stresses 
which  are  present  in  the  lead-through  seal  area  of  a microelectronic 
package  when  subjected  to  a thermal  shock. 

A' 

Subsequent  experiments  will  be  conducted  to  verify  the  predicted 
strengths  of  the  seal  and  how  it  is  influenced  by  the  manufacturing 
process,  and  establish  the  validity  of  the  assumptions  used  in  this  work. 

These  experiments  will  extend  the  present  work  to  include  microelectronic 
packages  with  flat  leads  and  test  fluids  other  than  water. 

The  results  of  this  program  will  be  used  by  RADC  to  develop  design 
guidelines,  for  use  in  conjunction  with  the  screening  procedures  of 
MIL-STD-883,  applicable  to  hybrid  and  integrated  microcircuit  packages  and 
electrical  feed  throughs  for  microwave  device  packages. 


PETER  F.  MANNO 
Project  Engineer 


CHAPTER  I 


INTRODUCTION 

In  order  to  ensure  the  reliability  of  microelectronic 
packages  which  are  used  by  the  Armed  Services,  these  pack- 
ages are  subjected  to  various  tests  which  simulate  the 
environments  which  the  packages  may  experience  in  service. 

The  test  methods  which  are  in  use  are  described  in  Military 
Standard  883B,  Test  Methods  and  Procedures  for  Microelectronics. 
The  present  study  concerns  a detailed  investigaton  of  the 
Thermal  Shock  Test  (Method  1011.2)  . 

As  described  in  MIL-STD-883B , the  purpose  of  the  thermal 
shock  test  is  to  determine  the  resistance  of  the  device  to 
sudden  exposure  to  extreme  changes  in  temperature  which  may 
be  encountered  in  equipment  which  is  operated  intermittently 
in  low  temperature  areas.  Physical  damage  may  be  experienced 
by  the  seals  of  a package  during  thermal  shock  as  a result  of 
thermal  stresses  which  accompany  the  extreme  temperature  dis- 
tribution which  must  exist  during  the  transient  heat  transfer 
process  associated  with  the  sudden  change  in  temperature. 

The  thermal  shock  test  apparatus  consists  of  temperature- 
controlled  baths  containing  fluids  appropriate  for  the  tem- 
perature extremes  which  are  desired.  The  least  severe  test 
level  consists  of  baths  of  ice  water  and  boiling  water  (Test 
Level  A) . The  most  severe  test  level  uses  baths  of  liquid 
nitrogen  and  silicon  oil  FC  70  or  UCON  100  with  temperatures 
of  -195°C  and  200°C  (Test  Level  I).  The  thermal  shock  test 
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! 


procedure  involves  bringing  the  device  to  equilibrium  tempera- 
ture in  the  low  temperature  liquid,  then  rapid1 v transferring 
it  to  the  high  temperature  liquid  and  allowing  it  to  come  to 
equilibrium,  and  then  continuing  this  back  and  forth  trans- 
ferring between  the  two  baths  for  a prescribed  number  of 
cycles,  usually  15. 

Thomas  [1]  conducted  an  experimental  study  of  packages 
subjected  to  thermal  shock  and  found  that  leakage  appeared  to 
be  directly  correlated  with  the  thermal  shock.  He  suggested 
that  thermal  stresses  in  the  vicinity  of  the  lead-through  seals 
may  be  large  enough  to  permit  leakage  during  the  thermal  shock 
which  might  go  undetected  during  subsequent  gross  leak  testing 
of  the  package.  Other  literature  [2]  indicates  that  stresses 
during  thermal  shock  testing  can  lead  to  fracture  of  the  glass 
seal . 

In  view  of  these  findings,  it  appears  desirable  to  have  a 
procedure  for  predicting  or  estimating  the  magnitude  of  the 
thermal  stresses  which  may  be  present  in  the  lead-through  glass 
seal  of  microelectronics  packages  which  are  subjected  to  thermal 
shock  testing.  The  predicted  stress  conditions  can  be  used  by 
the  individual  performing  the  thermal  shock  test  to  select  the 
test  level  which  is  consistent  with  the  level  of  stress  which 
is  desired  in  the  test  package.  An  analysis  of  thermal  stress 
for  conditions  of  thermal  shock  could  also  be  used  by  the  package 
designer  to  ensure  that  the  fabricated  package  will  perform 
satisfactorily  during  the  anticipated  environmental  thermal 
shock  conditions  that  it  would  be  subjected  to. 
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The  objective  of  the  present  study  is  to  develop  a mathe- 
matical model  for  predicting  the  transient  thermal  stresses 
that  occur  in  a typical  annular  lead-through  glass  seal  of  a 
microelectronics  package  during  a thermal  shock  test.  Since 
microelectronics  packages  typically  have  many  lead-through 
seals  and  since  the  package  geometry  varies  considerably  from 
one  case  to  the  next  it  is  necessary  to  develop  a simplified 
model  for  purposes  of  analysis.  The  approximate  model  which 
is  adopted  in  the  present  work  is  very  similar  to  the  one 
used  by  Borom  and  Giddings  [5]  which  focuses  attention  on  the 
stresses  in  the  vicinity  of  a typical  lead-through  seal.  The 
considerations  underlying  the  formulation  of  the  model  are 
as  follows. 

Consider  a typical  package  such  as  the  one  shown  in  Fig.  1. 
In  this  case  the  leads  all  pass  through  the  base  of  the  package. 
When  a package  such  as  this  is  subjected  to  a thermal  shock 
by  immersing  it  in  a hot  or  cold  liquid  bath,  the  lead  changes 
temperature  most  rapidly  due  to  its  small  mass  and  high  thermal 
conductivity.  The  base  responds  somewhat  less  rapidly  because 
of  its  larger  mass.  Both  the  lead  and  the  base  are  aenerally 
made  of  Kovar  and  therefore  have  a high  thermal  diffusivity. 
Because  of  this  and  because  of  the  nature  of  the  heat  transfer 
process  from  the  liquid  to  the  lead  and  the  base,  it  is  quite 
reasonable  to  assume  that  the  temperature  of  the  lead  and  the 
temperature  of  the  base  are  uniform  throughout.^"  On  the  other 

For  justification  we  can  consider  an  infinitely  large  Kovar 
plate  of  0.04  in.  thickness,  initially  at  32°F  and  exposed  to 
a 212°F  fluid  on  one  side  while  the  other  side  is  insulated. 
Using  the  non-dimensional  temperature  vs.  non-dimensional  time 
charts  provided  by  Schneider  [7] , one  finds  the  temperature- 
time histories  of  the  two  surfaces  of  this  plate  to  be 
essentially  identical. 
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hand,  the  annular  glass  seal  has  a much  lower  thermal  diffu- 
sivity  and  has  very  little  area  exposed  to  the  liquid  bath 
environment.  The  temperature  distribution  in  the  glass  seal 
can  be  approximately  modelled  as  having  a radial  variation 
from  the  temperature  of  the  lead  on  the  inner  annular  surface 
to  that  of  the  base  on  the  outer  annular  surface.  The  axial 
variation  in  temperature  is  neglected  in  view  of  the  small 
surface  area  in  contact  with  the  fluid  bath  in  comparison 
with  the  much  larger  areas  of  the  glass  seal  in  contact  with 
the  lead  and  the  base. 

It  is  anticipated  that  substantial  mechanical  stresses 
can  be  developed  in  the  lead,  and  the  base  and  the  glass  seal 
as  a result  of  the  time  varying  nonuniform  temperature  distribu- 
tion existing  in  the  package  during  thermal  shock.  Furthermore, 
it  is  suspected  that  the  state  of  thermal  stress  can  be  reasonably 
approximated  by  using  the  simplified  model  of  Borom  and  Giddings 
15]  which  focuses  attention  on  a typical  glass  seal.  The  model 
is  illustrated  by  Fig.  2 which  shows  a typical  glass  seal 
consisting  of  a lead,  the  glass  annulus,  and  an  annular  portion 
of  the  package  base.  The  outside  radius  of  the  base  annulus 

r is  selected  so  that  the  radial  stress  on  the  outer  annular 
o 

surface  could  be  assumed  to  be  zero.  As  illustrated  in  Fig.  2, 
the  radial  stress  at  the  free  edge  of  the  base  of  the  package 
would  be  expected  to  be  zero.  The  location  of  the  edge  of 
the  package  in  relation  to  the  lead  and  seal  therefore  defines 
a model  radius  r at  which  one  can  assume  zero  radial  stress. 


Free  edge  of 
package 


Glass  seal 


Lead 


Package  base 


Outer  periphery  of 
model  assumed  to  have 
zero  radial  stress 


Fig.  2.  Typical  glass  seal  consisting  of  lead,  annular  glass 
seal  and  part  of  package  base  in  the  form  of  an 
annulus . 
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7. 

The  model  therefore  consists  of  the  lead,  the  annular  glass 
seal  and  the  annular  segment  of  base  material  surrounding 
the  seal  as  illustrated  in  Fig.  3.  This  model  forms  the 
basis  of  the  thermal  stress  analysis  which  is  presented  in 
Chapter  II. 

In  order  to  carry  out  the  thermal  stress  analysis,  a tem- 
perature analysis  must  be  performed  in  order  to  determine: 
a)  the  uniform  temperature  of  the  lead,  b)  the  uniform  tempera- 
ture of  the  base,  and  c)  the  radial  temperature  distribution  of 
the  glass  seal,  all  as  functions  of  time.  The  model  which  is 
used  for  the  temperature  analysis  takes  into  account  the  entire 
amount  of  base  material  and  lead  material  in  the  package  and 
therefore  is  somewhat  different  than  the  model  depicted  in 
Fig.  3.  The  temperature  analysis  is  presented  in  Chapter  III. 

A computer  program  was  developed,  incorporating  the  tempera- 
ture distribution  model  and  the  thermal  stress  model.  It  can 
be  used  to  obtain  the  transient  thermal  stresses  occurring 
during  a thermal  shock  test.  This  analysis  was  applied  to  a 
commercially  used  microelectronic  package  in  which  the  leads 
and  the  base  were  made  of  Kovar  and  were  sealed  to  each  other 
with  7052  Corning  glass.  The  transient  thermal  stresses  developed 
in  the  seal  during  the  thermal  shock  test  level  A of  MIL-STD- 
883B  were  computed  and  presented  in  the  form  of  graphs.  The 
analysis  was  then  applied  to  packages  with  different  seal 
geometries  but  made  of  the  same  materials.  These  results  were 
non-dimensionalized  and  presented  in  the  form  of  graphs.  These 
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Pig.  3.  Model  of  the  lead-through  consisting  of  lead,  glass 
seal  and  a part  of  the  base  as  three  concentric 
cylinders . 
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graphs  can  be  used  directly  to  determine  the  thermal  stresses 
in  any  package  for  which  the  ratio  of  the  Young's  modulus  of 
the  metal  to  that  of  the  sealing  material  as  well  as  the  ratio 
of  the  coefficients  of  thermal  expansion,  the  ratio  of  the 
values  of  the  Poisson's  ratio  and  the  heat  transfer  property 
ratios  are  equal  to  the  corresponding  material  property  ratios 
hL  ( A./N)  P c 

— 7T- T-T2  = 96.459,  ^ — = 1.681,  E /Eo  = 2.439, 

k (L-L)  2 p c„  m se 

se  o se  se 

am/age  = 1-269,  vn/vse  = 1-364)  of  the  experimental  package 
used  in  this  investigation.  It  should  be  noted  that  since 
the  same  combination  of  materials  as  above  is  used  in  most 
microelectronic  packages,  the  present  study  can  be  of  general 
use  to  designers  and  screeners  of  packages. 
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CHAPTER  II 


THERMAL  STRESS  ANALYSIS 

The  temperature  analysis  of  the  model  lead- throucih  seal 
is  presented  in  Chapter  III.  The  present  chapter  treats  the 
analysis  of  stress  in  the  three  cylinder  model  which  was 
introduced  in  Chapter  I.  The  present  analysis  is  carried  out 
assuming  that  the  time-temperature  histories  in  the  cylindrical 
lead,  the  cylindrical  glass  seal  and  the  cylindrical  base  is 
known.  The  analysis  which  is  presented  below  follows  the  same 
basic  assumptions  as  those  made  by  Poritsky  [3]  in  his  analysis 
of  two  concentric  cylinders  subjected  to  uniform  temperature 
changes  (hence  stresses  arising  from  the  mismatch  of  thermal 
expansion  coefficients  of  the  glass  seal  and  metal  lead  and 
base  materials) . It  should  be  noted  that  Borom  and  Giddings  (5] 
extended  Poritsky' s [3]  analysis  to  include  three  concentric 
cylinders.  The  general  problem  of  the  thermal  stresses  in  a 
system  of  many  cylinders  with  nonuniform  temperature  distribution 
was  studied  by  Gatewood  [4],  The  analysis  presented  below  is 
essentially  the  same  as  that  presented  by  Gatewood  [4], 

The  thermal  stress  equations  are  obtained  by  writing  the 
general  equations  of  uncoupled  thermo-elasticity  for  a typical 
element  of  material  in  each  of  the  lead,  glass  seal  and  base 
cylinders.  The  solution  of  the  elasticity  equations  with  the 
appropriate  boundary  conditions  at  the  surfaces  cf  the  cylinders 
provides  the  desired  relations. 

The  equation  of  equilibrium  for  a typical  element  (cf. 

Fig.  4)  is 
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do  o o 

+ JLZ-±  = 0 
dr  r 


(1) 


Since  the  temperature  distribution  is  assumed  to  depend  only  on 
the  radius  r,  the  problem  is  ax i symmetric  and  only  one  equation 
of  equilibrium  is  necessary. 

In  any  isotropic,  homogeneous,  linearly  elastic  material, 
the  stresses  and  strains  accompanying  any  temperature  rise  6 
are  related  as  follows: 

c -a0=^-[o  - v(o , + o ) ] 

r E r <p  z 


e,  “ ae  = v t°A  - v(0  + a ) ] 

<p  E <p  r z 


e - a9  = — (o  - v(o  + o,)] 
z E z r p 


(2) 


where  e,  E,  v,  a are  respectively  the  strain,  the  Young's  modu- 
lus, the  Poisson's  ratio  and  the  coefficient  of  thermal  expan- 
sion. All  variables  except  E,  v and  a are  functions  of  radius 
and  time.  The  inversion  of  these  equations  yields 

E 


$ (l+vHl-iSr  [VGr  + (1'V)  V V£z  ' (1'V)  a6] 


0.  = 


1 1 1 — TTi — t — T [(l-v)e  + ve  , + ve  - (1+v)  ae] 

r (1+v)  (l-2v)  r <j>  z • 


a = 


(3a) 


(3b) 


The  strain-displacement  relations  for  an  axisymmetric  strain 
field  are 

du  ~ - — (4) 


£ = j 

r dr 


u 

" r 


where  u represents  the  displacement  in  the  radial  direction. 
Substituting  (4)  in  (3a)  and  (3b),  then  (3a)  and  (3b)  into  (1), 
one  obtains 


d^u  1^  du  u (1+v)  d6 
dr2  + r dr  " r7  = (1-v)  “ dr 


(5) 
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The  solution  of  (5)  is  of  the  form 


( 1+v)  1 

u = — r-  a — 

(1-v)  r 


0rdr  + C^r  + — 


(6) 


r . 

l 


where  and  C2  are  constants  of  integration.  r^  is  the 
inner  radius  of  the  particular  component.  Substituting  (6) 
in  (4)  and  (3)  yields 

C„ 


(1  + V)  rQ 
Lr  = (1-v)  a[6 


1 


0rdr]  + C,  - — 5- 
1 r2 


(1+v) 

£<P  " (1-v)  “ 


h \ 

1 r . 


r . 

l 


0rdr  + C + — ^ 
1 r+ 


(7) 


and 


Ea  1 

°r  ~ (1-v)  r2 


rr 


E ci  c2  VEez 

0rdr+  u^7  [u=2^r  " 7?1  + (i+vV(i-2v) 


a . = 


Ea 


(1-v)  r2 


0rdr  + 


r . 
1 


C C 

+ -£]  * 


vEe 


(1+v)  (l-2v) 


Ea0 


( 1+v) (l-2v)  (1-v) 


, (1-v) 

az  z (1+v) (l-2v) ' 


1 2VEC1 
Ea6  ‘l-v1  + (1+v) (l-2v) 


(8) 


Representing  the  integral 


rr 


r . 
1 


erdr  by  I(r^,r)  and  using  a sub- 


script j to  denote  each  of  the  three  concentric  cylinders  of  the 
model,  equations  (6)  and  (8)  can  be  written  for  each  component. 


(1+v.)  a.  C 

u . = — V -2-1  (r.,r)  +C..r+  — 1 

j (1-vJ  r 1 r 


3 


E.a. 

= _ _J_J k. 

'rj  (1-v.)  r2 
3 


q„J=-  <?  r.-r  77  1 (r  ,r)  + 


_ fr 

(1+v.)  1 (l-2v  . ) r 

3 3 


h 


(9) 


V .E  .E 

- ^ 2 


(1+v.) (l-2v  . ) 
3 3 


2v .E.C.  . . 

_]  + 3 3 L] _ E a 0 T 

'’zj  “j^z  1 (1+v  . ) ( l-2v  . ) (1+v  .)  (l-2v  .)  j j j (1-v.) 

3 3 3 3 3 


= E -e_  h 


(1-v.) 
_3_ 


(10) 


(11) 


— 1 

14 

where  j = 0,1,2  designate  respectively  the  base  made  of  Kovar, 
the  glass  seal  and  the  lead  made  of  Kovar.  (See  Fig.  3) 

Equations  (9),  (10),  (11)  have  two  unknown  constants  for 

each  cylinder.  For  the  three  cylinders,  the  total  number  of 
constants  is  six.  Therefore,  six  equations  are  needed  to 
determine  the  constants.  An  additional  equation  is  required 
to  evaluate  e z , which  is  assumed  to  be  constant.  The  equations 
derive  from  the  following  conditions: 

1)  The  radial  displacement  at  the  center  of  the  cylinders 
has  to  be  finite. 

2)  The  respective  radial  displacement  U£  and  u^  at  the  lead- 
glass  interface,  as  well  as  u^  and  Uq  at  the  glass-base  inter- 
face have  to  be  equal  to  each  other,  i.e.,  Uj  = u^  at  r = 
and  u^  = Uq  at  r = r^. 

3)  The  respective  radial  stresses  at  the  same  interfaces 

as  above  should  be  equal  to  each  other,  i.e.,o  = a atr=r 

r2  rl 

and  a -a  at  r = r. . 
ri  ro  1 

4)  The  outermost  boundary  should  be  free  of  radial  stress, 

i.e.,  = 0 at  r = Tq. 

o 

5)  The  total  force  in  the  axial  direction  should  be  zero, 

2 

i.e.  , y a . dA  = 0. 

i-0  J 21 

A 

These  conditions  yield  7 equations  and  7 unknowns.  The 
closed  form  solution  of  these  equations  would  be  too  complex 
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are  obtained  for  the  unknowns.  Hence  the  thermal  stresses  can 
be  obtained  at  each  increment  of  time.  The  results  of  this 
analysis  are  equivalent  to  B.E.  Gatewood's  analysis  of  thermal 
stresses  in  m + 1 concentric  cylinders  applied  to  the  case  of 
3 cylinders  [8]. 

The  analysis  presented  in  this  chapter  >_jnstitutes  the 
main  tool  for  the  evaluation  of  the  transient  thermal  stresses 
in  and  around  the  glass  seals. 


The  evaluation  of  the  thermal  stresses  requires  the 
knowledge  of  the  temperature  distribution  in  the  three  cylin- 
ders at  every  instant  of  time.  For  this  purpose  a temperature 
distribution  model  is  developed  which  takes  into  consideration 
the  heat  transfer  and  geometrical  characteristics  of  the  lead, 
the  glass  and  the  base  assembly. 

The  physical  characteristics  of  the  temperature  distribu- 
tion model  consist  of  a single  lead  surrounded  by  its  annular 
glass  seal  which  is  in  turn  surrounded  by  a "base"  which  has 
1/Nth  the  mass  and  1/Nth  the  outside  surface  area  of  the  total 
base.  N is  equal  to  the  number  of  leads  in  the  package.  The 
lead  and  the  base  are  considered  to  be  lumped  masses  and  as 
such  conduct  heat  instantaneously.  It  is  assumed  that  heat  is 
transferred  between  fluid  and  lead  and  fluid  and  base  but  not 
between  fluid  and  glass.  In  a heating  situation  (see  Fig.  5) 
heat  is  conducted  to  the  glass  from  the  lead  and  (it  will  turn 
out)  from  the  glass  to  the  base.  It  is  also  assumed  that  the 
temperature  distribution  in  the  glass  is  axisymmetric  at  every 
instant  of  time  and  varies  in  the  radial  direction  only.  The 
temperature  is  assumed  to  be  continuous  across  the  glass-metal  ' 
interfaces . 

The  assumption  of  a lumped  mass  for  the  lead  and  the  base 
can  be  justified  by  the  facts  that  the  metal  from  which  the  lead 
and  the  base  are  constructed  (usually  Kovar)  has  a high  thermal 
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Outside 

(Fluid) 


Fig.  5.  Diagram  representing  the  heat  transfer  to  the  seal 
when  the  seal  is  heating. 
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conductivity  and  the  masses  of  the  two  components  are  small. 

(See  also  footnote  on  p.  3)  The  reason  for  assuming  that  no 
heat  is  transferred  between  the  glass  and  the  fluid  is  that 
the  outside  surface  area  of  the  glass  is  very  small  in  relation 
to  the  outside  surface  areas  of  the  lead  and  the  base.  In  addi- 
tion, the  thermal  diffusivity  of  glass  is  about  one  tenth  that 
of  the  metal.  Therefore,  the  heat  transferred  between  the 
fluid  and  the  glass  is  small  compared  to  the  heat  transferred 
between  the  fluid  and  the  metal  and  then,  by  induction,  between 
the  metal  and  the  glass.  The  one  directional  heat  flow  from 
lead  to  glass  and  from  glass  to  base  in  Fig.  5 is  due  to  the 
lead  having  a larger  surface  area  per  unit  mass  than  the  base 

in  the  package  used  and  therefore  the  temperature  of  the  lead 
increases  more  rapidly  than  the  temperature  of  the  base.  Were 
the  geometrical  and/or  physical  characteristics  of  the  lead  and 
the  base  different  the  glass  could  receive  heat  from  the  base 
and  the  lead  or  from  the  base  only.  The  reason  for  temperature 
variation  only  in  the  radial  direction  is  that  no  heat  is  trans- 
ferred from  the  fluid  to  the  glass.  Also,  since  the  lead  and 
the  base  each  are  assumed  to  have  a uniform  temperature  distri- 
bution the  only  temperature  gradient  exists  in  the  radial  direction. 


The  temperature  distribution  model  is  used  to  determine  the 
transient  temperature  distribution  in  the  glass.  The  transfer 
of  heat  is  governed  by  the  partial  differential  equation. 


+ i|I) 

r 3r 


3T 

3t 


(12) 
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where  a = k /p  c is  the  thermal  diffusivity  and  k , p , c 
g 9 g 9 9 9 

are  respectively  the  thermal  condictivity , the  density  and  the 
specific  heat  of  glass.  t and  r indicate  respectively  the 
time  and  the  radial  coordinate.  T is  the  temperature  in  the 
glass,  which  is  a function  of  the  radius  and  time. 

The  initial  condition  is: 

T(r,0)  = Tq,  an  initial  uniform  temperature  for  the 
lead-base-glass  assembly  corresponding  to  the  beginning 
of  the  thermal  shock. 


k 

g ar1 


— — ^ (t  - t|  ) - p c !I| 

A N °°  'r=r,  pk  k A.N  9t! 
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where  Tm  is  the  temperature  of  the  fluid  and  h is  an  effective 
heat  transfer  coefficient  governing  the  rate  of  heat  transfer 
between  the  lead  or  the  base  and  the  fluid.  h is  determined 
through  the  experimentally  measured  temperature-time  histories 
for  the  lead  and  the  base  as  described  below.  The  as  yet 
undefined  variables  in  equations  (14)  and  (15)  are  defined 
as  follows: 


Pk  = the  density  of  the  metal  (usually  Kovar) 

ck  = the  specific  heat  of  the  metal 

A = the  outside  surface  area  of  the  lead  in  contact 

with  the  fluid 


j 


m — "UJI 


= the  outside  surface  area  of  the  base  and  lid 

A^  = the  surface  area  of  the  glass-base  interface 

A£  = the  surface  area  of  the  lead-glass  interface 

V = the  total  volume  of  the  lead 

a 

V,  = the  total  volume  of  the  base  and  lid  ^ 
b 

N = the  total  number  of  leads. 

The  differential  equation  (12)  and  boundary  conditions 
(13,14,15)  are  solved  using  a finite  element  analysis.  For 
this  purpose  the  glass  cylinder  is  divided  into  5 radial  elements 
of  equal  radial  width.  The  physical  properties  of  glass  are 
assumed  to  be  uniform.  Details  of  the  analysis  are  described 
in  Appendix  A. 

The  solution  of  equation  (12)  yields  the  transient  tem- 
perature distribution  in  the  glass  as  well  as  the  temperature- 
time histories  of  the  lead  and  the  base. 


The  microelectronic  package  used  in  this  investigation  consists 
of  a base  and  a lid  as  shown  in  Fig.  8.  Other  packages  have  the 
base  surrounded  with  vertical  walls  with  a lid  in  the  form  of  a 
plate.  For  these  types  of  packages  the  walls  should  be  included 
in  the  analysis.  More  generally,  a^  and  Vb  refer  to  all  the 
parts  of  the  package  except  the  glass  seals  and  the  leads. 


21. 


CHAPTER  IV 
EXPERIMENTS 

4.1  Procedure 

The  approximate  radial  temperature  distribution  through- 
out the  3 cylinder  model  can  be  determined  analytically  if 
the  temperature  of  the  lead  and  base  cylinders  is  known  as 
a function  of  time.  These  temperature  histories  can  be 
determined  if  the  effective  heat  transfer  coefficient  for  the 
lead  and  base  are  known.  In  the  present  work,  the  temperature- 
time histories  of  the  lead  and  the  base  of  the  package  were 
measured  experimentally.  This  data  was  then  used  to  evaluate 
the  effective  heat  transfer  coefficient  h which  represents  the 
rate  of  heat  per  unit  area  transferred  to  the  package  during 
given  experimental  conditions. 

The  standard  procedure  for  performing  thermal  shock  tests 
as  given  in  MIL-STD-883B  [6]  requires  that 

(a)  The  package  be  immersed  in  the  low  temperature  bath 
for  a minimum  of  5 minutes. 

(b)  The  package  be  transferred  from  the  low  temperature 
to  the  high  temperature  bath  and  left  there  for  5 minutes. 

(c)  The  package  be  transferred  from  the  high  temperature 
to  the  low  temperature  bath  and  left  there  for  5 minutes.  The 
remaining  stems  are  the  repetition  of  b and  c at  least  fourteen 


times . 
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The  experiments  which  were  performed  consisted  of  recording 


the  temperature  of  the  lead  and  the  base  as  a function  of  time 


during  and  subsequent  to  immersion  of  the  package  in  the  thermal 
shock  fluid  bath.  Experiments  were  carried  out  for  level  A 
conditions  of  ice  water  and  boiling  water  baths. 

The  packages  used  during  the  experiments  /ere  Tekform  packages 
identified  with  part  number  20133,  style  A in  the  manufacturer's 
catalog.  A schematic  representation  of  the  packages  is  given 
in  Fig.  6. 

The  temperature  was  measured  using  Chromel-Constantan 
thermocouples.  The  hot  junction  of  the  thermocouples  was  formed 
by  twisting  the  two  wires  together.  The  thermocouple  wires  were 
of  .005  in.  diameter  and  about  2 ft  long.  They  were  soldered 
to  the  lead  and  the  base  with  rosin  fluxed  66  Sn  - 34  Pb  solder. 

(See  Figures  7a  and  7b.)  The  thermocouples  were  connected  to 
a multichannel  Honeywell  Visicorder  (Model  1108)  instrumented 
with  M40  350A  type  Honeywell  galvanometers.  The  Visicorder  was 
used  to  simultaneously  record  the  output  e.m.f.  of  each  of  the 
thermocouples.  The  experimental  set  up  is  shown  in  Fig.  8. 

No  controlled  reference  cold  junction  was  used  during  the  experiments. 
Instead,  the  recorder  was  set  to  an  initial  position  corresponding 
to  the  lower  temperature  bath.  The  calibration  of  the  recorder 
and  thermocouples  was  then  used  to  determine  the  temperature  corres- 
ponding to  any  voltage  output.  Also,  the  duration  of  the  heating  or 
cooling  process  was  very  short  (7-8  sec).  Therefore  variation  of 
temperature  of  the  room  had  a negligible  effect  on  the  measurements. 


Fig.  6.  Schematic  representation  of  the  Tekform 
package  used  during  experiments. 


Photograph  showing  the  thermocouple 
on  the  base  and  the  thermocouple  on 
lead.  Note  the  solder  on  the  four 
sides  of  the  base. 


Fig.  8.  Photograph  showing  the  set-up  for  the 
experiment.  From  left  to  right: 
Visicorder,  ice  water  container; 
boiling  water  on  heater  with  the 
package  in  it. 
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Using  this  set  up,  the  package  was  immersed  in  ice  water 
( 32°F)  and  left  there  for  about  5 minutes.  Then,  the  recorder 
was  activated  as  the  package  was  transferred  to  boiling  water 
(212°F) . After  a few  seconds  the  recorder  was  stopped.  By 
this  time  the  package  had  reached  about  90  percent  of  its 
equilibrium  temperature.  After  complete  thermal  equilibrium, 
the  package  was  immersed  back  into  ice  water  and  the  cooling 
process  recorded. 

The  main  problems  encountered  with  performing  the  experi- 
ments were  related  to  leakage  of  the  package.  The  lid  was 
sealed  to  the  base  initially  using  epoxy  glue.  After  two 
repetitions  of  the  cycle  the  lid  seal  would  weaken  causing 
water  to  leak  into  the  package,  which  in  turn,  caused  distur- 
bances in  the  recording  of  the  thermocouple  e.m.f. 

This  difficulty  was  overcome  by  using  solder  instead  of 
glue  to  seal  the  lid  to  the  base.  The  method  by  which  this 
was  accomplished  can  be  briefly  described  as  follows: 

The  lid  is  placed  on  a hot  plate  and  brought  to  a tem- 
perature of  approximately  175°C.  The  base  is  then  placed  on 
the  lid  and  brought  to  the  same  temperature.  A heat  sink 
consisting  of  copper  blocks  is  placed  on  the  base,  preventing 
the  solder  attaching  the  thermocouple  wires  from  melting.  A 
fine  bead  of  Indalloy  solder  (the  melting  temperature  of  which 
is  lower  than  that  of  66Sn  - 34Pb  solder)  is  put  on  the 
boundaries  of  the  base  and  lid.  Finally  a soldering  iron  is 
run  gently  along  the  sides,  causing  the  Indalloy  solder  to  melt. 
The  whole  assembly  is  then  left  to  cool. 
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The  package  thus  sealed  was  tested  for  hermeticity  using  the 
following  procedure.  The  package  was  subjected  to  a pressure  of 
10  psi  for  30  minutes  in  a pressure  bomb.  (The  applied  pressure 
of  10  psi  was  calculated  not  to  result  in  any  bond  failure  [9]). 
The  pressurized  package  was  then  immersed  in  warm  water  and 
checked  for  air  bubbles. 


4.2  Evaluation  of  the  Effective  Heat  Transfer  Coefficient 

The  experimental  temperature-time  histories  of  the  lead 
and  the  base,  recorded  using  the  Visicorder  are  shown  in 
Fig.  9a,  b,  c,  d,  e,  f,  g and  h\  These  curves  were  used  to 
evaluate  the  effective  heat  transfer  coefficient  h. 

In  the  thermal  model,  presented  in  the  preceding  chapter, 
h determines  the  rate  of  heat  transfer  to  the  lead  and  the 
base.  h appears  in  the  boundary  conditions  of  the  partial 
differential  equation  governing  the  heat  transfer  to  the  glass. 

When  this  differential  equation  is  solved,  the  temperature-time 
histories  corresponding  to  the  lead  and  the  base  are  obtained. 

The  experimental  value  of  h is  found  in  such  a way  as  to  ensure 
that  the  mathematically  determined  temperature-time  histories 
provided  best  fit  of  the  experimental  temperature-time  histories. 
The  best  fit  is  realized  by  solving  the  differential  equation 
for  different  values  of  h , until  the  h value  yielding  the  smallest 
least  square  error  between  the  values  corresponding  to  the 
first  2.5  seconds  of  the  mathematical  and  experimental  temp- 
erature-time curves  is  obtained.  The  least  square  error  is 

n 2 

defined  by  Er  = . £ (y^  - Ye^)  » where  y^  and  yg^  are,  respect- 

ively corresponding  points  on  the  theoretical  and  experimental 
time-temperature  histories.  A typical  example  illustrating 


1 The  code  on  these  figures  is  explained  on  pace  41 
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Fig.  9e  Experimental  temperature-time  histories  for 

the  lead  and  the  base  obtained  from  Experiment 
3 when  the  package  at  32 °F  is  immersed  into  a 
212 °F  bath. 
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the  correspondence  between  the  theoretical  and  experimental 
curves  is  shown  in  Fig.  10. 

The  least  square  fit  over  the  first  2.5  seconds  of 
the  experimental  and  theoretical  time-temperature  histories 
was  done  for  the  following  reasons.  When  the  experiments 
were  performed,  after  immersing  the  package  from  the  cold 
bath  into  the  hot  bath,  the  cold  package  caused  the  water  to 
cool  to  a certain  extent}  It  is  felt  that  this  behavior 
would  not  affect  the  first  2.5  seconds  of  the  experiment. 

It  will  also  be  shown  in  Chapter  V that  the  maximum  transient 
stresses  occur  during  the  first  interval  of  the  thermal  shock 
test . 

As  previously  stated,  the  experiments  were  performed 
using  thermocouples  with  thermocouple  wires  soldered  to  the 
lead  and  the  base.  The  additional  mass  consisting  of  the  thermo- 
couple wires  and  solder  is  negligible  compared  to  the  mass  of 
the  base.  However,  when  the  lead  is  considered  the  presence  of 
the  thermocouple  wires  and  the  solder  becomes  important  since  these 
masses  are  comparable  in  magnitude  to  that  of  the  lead.  Taking 
into  account  the  mass  of  the  thermocouple  and  the  solder  the  boun- 
dary condition  (14)  of  the  partial  differential  equation  (12) 
becomes 
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1 The  theoretical  and  experimental  curves  in  Fig.  10  diverge 
because  the  theoretical  curve  reaches  212°F  while  the 
experimental  curve  remains  at  a lower  value  due  to  the 
cooling  of  the  bath. 
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where  the  additional  subscripts  s,  ch,  co  represent  respectively 
the  solder,  the  chromel  thermocouple  wire  and  the  constantan 
thermocouple  wire.  The  solder  bulk  is  assumed  to  be  a sphere  of 
0.073  in.  diameter,  an  approximate  value  obtained  from  the  actual 
solder  on  the  package.  For  calculation  purposes  1 in.  of  each 
of  the  thermocouple  wires  is  assumed  to  affect  the  heat  transfer 
process.  The  length  of  thermocouple  wire  to  be  included  is 
completely  unknown  and  only  an  appropriate  estimate  can  be  made. 
The  length  of  each  wire  inside  the  package  is  approximately  1 in. 
Therefore  this  amount  was  assumed  to  affect  the  heat  transfer 
process.  Thus  the  lead  and  the  attached  components  can  be 
visualized  as  shown  in  Fig.  11. 

It  should  be  noted  that  these  additional  terms  do  not 
affect  the  first  term  of  the  right-hand  side  of  equation  (14b). 
This  is  because  the  thermocouple  wires  and  solder  are  all  in- 
side the  package  and  do  not  have  a direct  contact  with  the 
fluid . 

Thus  when  h is  evaluated  using  the  experimental  curves, 
equation  (14)  is  replaced  by  equation  (14b)  as  the  boundary 
condition  in  the  mathematical  model. 

The  numerical  values  which  apply  to  the  experimental  pack- 
age and  used  in  the  calculation  of  the  coefficients  in  equation 
(12),  (14),  and  (14b)  of  the  temperature  distribution  model 


are : 
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k = 1.366  x 10"5  Btu/ sec-in- °F 

g 

p = 8.2  x 10~2  lb/in3 

g 

Pk  = 0.302  lb/in3 

Ps  = 0.304  lb/in3 

p . = 0.315  lb/in3 

ch 

p = 0.322  lb/in3 
co 

C = 0.23  Btu/lb- °F 

g 

ck  = 0.105  Btu/lb- °F 

c = 0.051  Btu/lb- °F 
s 

C . = 0.107  Btu/lb-°F 

ch 

c = 0.094  Btu/lb- °F 
co 

L = 0 . 205  in 
o 

L = 0.295  in 

L = 1 in 
v? 


d = 0 . 018  in 


0.073  in 

d = 0 . 005  in 
co 

0.15434605  in3 
4.2874  in2 
0.000075068  in3 
0.0115925  in2 


N = 48 
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4.3  Expe rime ntal  Results 

The  effective  heat  transfer  coefficient  was  determined 
from  results  of  experiments  carried  out  with  three  Tekform 
packages  which  were  submitted  to  thermal  shock  tests.  The 
packages  were  prepared  and  tested  using  the  experimental 
technique  described  in  the  preceding  section.  The  first  two 
packages  were  sealed  with  glue  and  the  third  with  solder. 

The  effective  heat  transfer  coefficients  which  were  evaluated 
using  the  experimentally  determined  temperature-time  histories 
are  tabulated  in  Table  1. 

Each  row  in  Table  1 corresponds  to  a single  experiment  (i.e.. 
Experiment  1 indicates  that  the  first  package  was  submitted  to 
one  thermal  shock).  The  experiment  numbering  code  is  as  follows: 

H denotes  a heating  experiment  (package  at  32°F  transferred  to  a 
212°F  bath) ; C denotes  a cooling  experiment  (package  at  212°F 
transferred  to  a 32°F  bath).  1,  2,  and  3 denote  the  package 

t 

numbers.  (a)  and  (b)  signify  the  first  and  2nd  experiment  of 
the  designated  type  if  these  were  more  than  one.  E.g,  Experiment 
C3 (b)  is  the  second  cooling  experiment  performed  with  package  3. 

t 

The  columns  in  Table  1 refer  to  the  location  of  the  ther- 
mocouple on  the  package.  The  first  column  represents  the 
effective  heat  transfer  coefficient  h calculated  from  a tem- 
perature-time history  obtained  from  a thermocouple  attached 
to  the  lead.  Columns  (2)  and  (3)  are  the  effective  heat 
transfer  coefficients  h for  the  package  as  determined  by 
temperature  measurement  at  the  center  of  the  base  and  at  the 
corner  of  the  base  respectively.  Values  for  both  columns  (2) 
and  (3)  represent  the  h for  the  base  since  in  the  temperature 
distribution  model,  the  base  was  assumed  to  be  one  lumped  mass 
having  a uniform  temperature  at  each  instant  of  time. 


! 
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It  should  be  noted  that  in  Table  1,  certain  values  of  h 
(e.g.,  Experiment  2,  "cooling")  are  not  given.  This  is  because 
water  leakage  through  the  lid  seal  of  the  package  during  these 
experiments  rendered  the  thermocouple  voltage  of  the  recorder 
unreadable.  Thus,  the  value  of  h for  these  cases  could  not  be 
determined . 

The  effective  heat  transfer  coefficient  of  an  object 
depends  on  the  convection  process  and  therefore  also  depends 
on  its  geometry,  its  position  in  the  surrounding  environment 
and  the  process  by  which  heat  is  transferred  to  or  from  it. 

It  is  conceivable  therefore  that  the  true  value  of  h for  the 
lead  differs  from  that  of  the  base  since  the  conditions 
mentioned  above  are  different  for  the  lead  and  the  base.  This 
difference  can  be  explained  by  use  of  Fig.  12(a)  and  12(b). 

Fig.  12(a)  represents  a cold  lead  situated  in  a hot  fluid. 

The  entire  surface  of  the  cold  lead  is  continually  in  contact 
with  new  particles  of  hot  fluid.  Fig.  12(b)  shows  the  cold 
base  standing  in  the  hot  fluid.  A point  A on  the  bottom  part  of 
the  base  is  in  contact  with  the  hot  fluid  which  transfers 
heat  to  this  part  and  moving  towards  B becomes  colder.  So  B 
is  in  contact  with  cooler  particles  of  fluid.  This  heating 
process  would  cause  the  lead  to  heat  faster  than  the  base  which 
would  result  in  a higher  effective  heat  transfer  coefficient 


for  the  lead  than  for  that  for  the  base. 


45. 


However,  as  explained  before,  the  assumption  of  a certain 
mass  of  solder  and  an  estimated  length  for  the  thermocouple 
wires  results  in  considerable  uncertainty  in  the  calculation  for 
the  effective  heat  transfer  coefficient  for  the  lead.  Fur- 
thermore, the  experimentally  measured  values  of  h shown  in 

2 

Table  1 vary  from  70  to  375  Btu/ft  -hr-°F,  a rather  large 
range.  Due  to  this  large  range  in  experimentally  determined 
values  coupled  with  the  uncertainty  in  the  calculation  for  the 
lead  it  is  difficult  to  justify  the  assumption  that  the  heat 
transfer  coefficient  for  the  lead  is  different  from  that  of  the 
base.  The  evaluation  of  R for  the  base  is  considered  to  be  more 
accurate  , therefore  the  average  value  of  R determined  for  the  base 
is  taken  to  be  the  most  reliable  estimate  of  the  overall  effective 
heat  transfer  coefficient  for  both  the  lead  and  the  base. 

The  transient  temperature  distribution  in  the  glass  seal 
which  consists  of  the  solution  of  Eq.  (12)  and  boundary  condi- 
tions (14)  and  (15)  can  be  obtained  when  h is  known.  The  tran- 
sient temperature  distribution  in  the  glass  seal  of  the  experi- 
mental package,  thus  obtained,  is  shown  in  Fig.  13  and  14.  (It 
should  be  noted  that  these  temperature  distributions  correspond 
to  the  package  without  any  attachment  (thermocouple  wires  and 
solder)  in  it) . 
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Radial  temperature  distribution  in  the  glass  annulus  at 
different  increments  of  time  when  the  seal  is  cooled 
from  2 12 °F  to  32°F. 
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CHAPTER  V 
RESULTS 

The  results  of  the  thermal  stress  calculations  for  the 
Tekform  package  which  was  used  in  the  experimental  portion  of 
the  study  regarding  the  measured  effective  heat  transfer  coef- 
ficient are  presented  in  this  chapter.  The  results  are  presented 
in  such  a way  as  to  apply  to  other  packages  providing  they  are 
characterized  by  having  certain  ratios  of  elastic  modulus,  co- 
efficient of  thermal  expansion  and  Poisson's  ratios  for  the  metal 
and  glass  materials  of  the  package  as  described  in  detail. 

The  calculation  of  the  stresses  was  performed  by  means  of 
an  APL  program  listed  in  Appendix  B.  The  program  accepts 
as  input  the  material  properties  (Young's  modulus  (psi) , coeffi- 
cients of  linear  expansion  (in/in/°F),  Poisson’s  ratio),  the  three 
radii  corresponding  to  the  three  cylinders  (in),  the  respective 
densities  of  Kovar  and  glass  (lb/in2) , the  respective  specific 
heats  of  Kovar  and  glass  (Btu/lb-°F) , the  thermal  conductivity 
of  glass  (Btu/in-sec-°F) , the  outside  surface  area  of  the  lead 
in  contact  with  the  fluid  (in2),  the  outside  surface  area  of 
the  base  and  lid  (in2) , the  volume  of  one  lead  (in^) , the  total 

1 n # 

volume  of  the  base  and  lid  (in-5),  the  initial  temperature  of 
the  package  before  the  thermal  shock  test  (°F) , the  final  tempera- 
ture the  package  reaches  after  the  test  (°F) , the  effective  heat 
transfer  coefficient  h (Btu/hr-f t^-°F) , the  radius  where  the 
stresses  are  to  be  evaluated  (in),  the  total  number  of  leads, 

1 See  footnote  on  page  20. 
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the  time  increment  (sec),  and  the  total  number  of  time  incre- 
ments. The  output  consists  of  the  transient  absolute  and  non- 
dimensional  radial  and  tangential  stresses  developed  at  the 
given  radius.  A typical  example  illustrating  the  use  of  the 
program  and  the  output  is  given  in  Appendix  C. 

The  results  are  presented  in  terms  of  the  following  non- 
dimensional  parameters: 
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radial  parameter 
temperature  ratio 
time  parameter 

effective  heat  transfer  coeffi- 
Young's  modulus  ratio 


a = non-dimensional  thermal  expansion  coefficient 

ratio 

a*  = non-dimensional  stress 

0 = difference  between  the  high  and  low  extreme 

temperatures  of  a thermal  shock  test 

k = subscript  representing  Kovar 

g = subscript  representing  glass 

rr  = 2.05  sec.  which  is  the  time  necessary  to  bring 

to  90%  of  its  equilibrium  temperature  a cube  of 
Kovar  with  a mass  of  1 lb  when  h = 210/Btu/hr-f t2- 
°F  while  heating  from  32°F  to  212°F. 

When  the  thermal  stress  analysis  is  applied  to  the  experi- 
mental package  the  results  obtained  are  presented  in  non-dimensional 
form.  It  should  be  noted  that  Fig.  15a  through  Fig.  18  are  the 
illustration  of  the  results  for  the  specific  package  used  in  this 
investigation.  Graphs  shown  in  Fig.  19  through  25  apply  to 
packages  with  circular  leads  and  with  the  lead  and  the  base  made 
of  the  same  material.  Also  to  be  able  to  use  these  graphs,  the 
ratio  of  the  Young's  modulus  of  the  metal  to  that  of  the  sealing 
material  for  the  package  under  consideration,  as  well  as  the 
ratio  of  the  coefficients  of  thermal  expansion  and  the  ratio  of 
the  values  of  the  Poisson's  ratio  must  be  equal  to  the  correspond- 
ing material  property  ratio  of  the  experimental  package  (viz., 

E./E  =2.439;  a, /a  = 1.269;  v, /v  = 1.364  . For  Fig.  19  and  20  which 
K g K g k g 

represent  the  maximum  transient  stresses,  the  conditions  that 
hL(Ab/N)/kg(L-L0) 2 = 96.459  and  Pkck/pgcg  = 


1.681  must  be  satisfied. 
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When  the  Tekform  package  was  taken  from  32°F  to  212°F  the 
resulting  transient  non-dimensional  radial  and  tangential  stresses 
were  computed  and  plotted  in  Figs.  15a,  15b,  16a,  16b.  The  di- 
mensional values  of  the  stresses  can  be  obtained  for  the  stresses 
in  these  graphs  by  using  Eq.  22  which  can  be  written  as  a = 
o*EgZ  The  dashed  lines  in  each  of  these  figures  represents 

the  stress  distribution  in  the  3-cylinder  model.  The  3 regions 
are  identified  on  the  graphs  as  the  lead,  the  glass  and  the  base. 
Each  line  shows  the  stress  distribution  at  a different  time. 

The  direction  of  the  change  in  the  stress  distribution  as  time 
increases  is  shown  on  each  of  these  figures. 

It  is  assumed  in  these  graphs  that  at  t = 0 the  entire  pack- 
age is  stress-free.  Fig.  15a  shows  that  as  time  increases  the 
radial  stresses  become  compressive  for  all  values  of  radius  and 


the  magnitude  of  the  compressive  stress  at  each  radius  increases 
until  t =0.9  sec.  (At  this  time  the  maximum  non-dimensional 
compressive  radial  stress  at  the  lead-glass  interface  is 
-0.617  which  corresponds  to  -2333  psi)  Fig.  15b  represents 
the  continuation  of  the  process.  The  graph  shows  that  the 
magnitude  of  the  compressive  radial  stress  at  each  radial 
position  decreases  as  time  increases.  When  finally  thermal 

equilibrium  is  reached  at  t - 12  sec,  the  radial  stress 
achieves  a relatively  small  tensile  value  at  the  lead-qlass 

interface  and  at  the  glass-base  interface.  Inspection  of 
Figs.  15a,  and  15b  shows  that  the  largest  radial  stress  magni- 
tude occurs  at  the  lead-glass  interface  soon  after  immersion 
in  the  hot  fluid  bath.  Investigation  of  several  examples 
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Fig.  15b  Non-dimensional  radial  stress  at  each  radius 
as  a function  of  time. 
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having  different  values  of  rQ,  r^,  , showed  that  this 

general  behavior  regarding  the  location  of  the  highest  magni- 
tude of  stress  at  the  lead-glass  interface  was  invariably 
obtained.  In  subsequent  analysis  it  has  been  assumed  that 

the  most  severe  radial  stress  always  occurs  at  the  lead-glass 
interface.  Of  course,  the  most  severe  stressing  from  the  fail- 
ure point-of-view  would  be  expected  during  immersion  in  the 
ice  bath  which  would  result  in  large  tensile  radial  stress  at 
the  lead-glass  interface. 

The  same  process  is  shown  in  Fig.  16a  and  16b  where  the 
tangential  stress  distribution  in  the  heating  process  is  plot- 
ted for  different  times.  Fig.  16a  shows  that  a small  portion 
of  the  glass  near  the  glass-base  interface  experiences  compres- 
sive stress.  However,  for  most  radial  positions  the  tangential 
stress  is  tensile  and  the  magnitude  of  the  stresses  at  each 
radial  position  becomes  larger  as  time  increases.  Fig.  16b 
shows  that  the  tangential  stress  at  each  radius  decreases  at 
the  right  hand  side  of  the  glass-base  interface  and  increases 
at  the  left  hand  side  of  the  glass-base  interface  until 
t 12  sec  when  the  temperature  of  the  entire  package  is 
equal  to  the  temperature  of  the  fluid.  This  process  can  be 
seen  better  in  Fig.  17  and  18  where  the  non-dimensional  radial 
and  tangential  stresses  at  the  lead-glass  and  glass-base  inter- 
face are  plotted  as  a function  of  the  non-dimensional  time. 

Fig.  17  shows  that  the  radial  stress  at  the  lead-glass  interface 
decreases,  reaches  a minimum  and  subsequently  increases  until  the 
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Fig.  16a  Non-dimensional  tangential  stress  at  each  radius 
as  a function  of  time. 


-dimensional  tangential  stress  at  the  glass-base  and  lead-glass 
erf ace  as  a function  of  non-dimensional  time. 
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reached.  Thenceforth,  the  radial  stress  increases  until  its 

equilibrium  value  is  achieved.  These  maximum  and  minimum 

stresses  occur  when  the  temperature  difference  across  the  3 

cylinders  becomes  the  largest.  As  the  heating  process  continues 

this  gradient  decreases  causing  the  radial  stresses  to  diminish. 

When  the  steady-state  temperature  is  reached  the  radial  stresses 

come  to  a steady-state  value.  This  residual  stress  is  caused 

by  the  mis-match  in  thermal  expansion  coefficients  of  Kovar  and 

glass  for  this  temperature  range  (a^  = 3.25  x 10  6 in/in- °F, 

a = 2.56  x 10  6 in/in- °F  for  32°F  < T < 572°F)  . P.ecall  that  it 

y 

was  assumed  that  the  stresses  are  zero  at  the  beginning  of  the 
heating  period. 


i 


The  transient  tangential  stress  at  the  lead-glass  inter- 
face is  shown  in  Fig.  18.  This  stress  passes  also  through 
a maximum,  then  decreases  to  come  to  steady-state.  The  absolute 
value  of  the  maximum  tangential  stress  in  this  case  is  1685  psi. 

On  the  other  hand,  the  tangential  stress  at  the  glass-base 
interface  decreases  slightly  to  reach  minimum  value  and  then 
increases  until  it  reaches  the  steady-state  value.  These  results 
illustrate  the  stress  behavior  of  the  seal  when  it  is  submitted 
to  the  heating  cycle  of  the  thermal  shock  test.  In  the  case  of  a 
cooling  cycle,  the  numerical  results  for  all  of  the  stresses 
would  be  the  same  as  those  illustrated  for  a heating  cycle,  how- 
ever, the  stresses  would  be  of  opposite  sign. 

In  order  to  assess  the  importance  of  knowing  the  true  value 
of  the  overall  effective  heat  transfer  coefficient  h for  the 
package  insofar  as  the  calculated  maximum  stress  is  concerned, 

1 
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calculations  were  made  for  a range  of  values  of  h consistent 
with  the  probable  range  of  uncertainty  as  determined  by  experi- 
ment. As  described  in  section  4.3,  the  average  value  of  h was 

found  to  be  210  Btu/f t^-hr-°F,  however,  the  experimental  values 
ranged  from  a low  of  70  Btu/f t2-hr-°F  to  a high  of  375  Btu/ft^- 
hr-°F.  The  influence  of  changing  the  value  of  h in  the  calcula- 
tion for  the  stresses  is  illustrated  in  Table  2 which  records 
the  maximum  value  of  non-dimensional  radial  stress  (at  the  lead- 
glass  interface)  for  varying  values  of  h ranging  from  130-290  Btu/ 
ft  -hr-°F.  As  shown  in  the  table,  the  maximum  difference  from 
the  stress  calculated  from  210  Btu/f t^-hr-°F  was  about  10  percent. 
This  is  not  judged  to  be  a large  difference,  and  it  is  felt  that 
some  uncertainty  in  the  true  value  of  h does  not  result  in  a very 
large  uncertainty  in  the  maximum  value  of  the  calculated  radial 
stress  at  the  lead-glass  interface.  On  the  other  hand,  calcula- 
tions have  shown  that  if  h is  assumed  to  be  ten  times  larger  than 
its  average  value  (210  Btu/f t2-hr-°F)  the  maximum  radial  stress 
becomes  about  37  percent  larger  and  if  h is  ten  times  smaller 
than  its  average  value  the  maximum  radial  stress  is  about  48  per- 
cent smaller  than  the  maximum  radial  stress  corresponding  to  the 
average  value  of  h.  This  represents  an  important  variation.  How- 
ever an  effective  heat  transfer  coefficient  which  is  that  much 
different  than  the  average  value  would  correspond  to  a different 
heat  transfer  environment.  In  fact,  the  heating  and  cooling  pro- 
cedure used  in  thermal  cycling  tests  (as  opposed  to  thermal  shock 
tests)  is  characterized  by  a very  slow  heat  transfer  process  for 
which  h would  be  expected  to  be  very  small  compared  to  the  h 


h (Btu/f t2 
'"-^hr-°F) 
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obtained  in  the  thermal  shock  experiments.  Therefore,  it  would 
be  reasonable  to  conclude  that  the  small  stresses  evaluated 
using  an  effective  heat  transfer  coefficient  10  times  smaller 
than  the  average  experimentally  obtained  h could  be  indicative 
of  the  stresses  occurring  during  a temperature  cycling  test. 

The  analysis  concerning  the  uncertainty  of  h was  based  on 
the  assumption  that  the  lead  and  the  base  have  an  equal  effec- 
tive heat  transfer  coefficient.  In  Chapter  IV,  however,  it 
was  shown  that  there  is  some  justification  for  believing  that 
the  lead  and  the  base  could  have  different  values  of  h.  The 
effect  of  such  variation  was  investigated  by  assuming  values  of 
h for  the  lead  and  the  base  differing  by  20,  40  and  60  Btu/ 
ft2-hr-°F  and  calculating  the  thermal  stresses  in  each  case. 

The  values  of  h given  to  the  base  and  the  lead  were  always  cen- 
tered around  the  average  value  of  h obtained  experimentally. 

The  results  are  tabulated  in  Table  3.  The  percent  change  between 
the  stresses  evaluated  using  different  h's  (for  the  lead  and 
the  base)\  and  the  stresses  obtained  using  the  average  h for 
both  of  them  is  shown  to  have  a maximum  value  of  15.5  percent. 
This  is  not  judged  to  be  a large  difference.  This  result 
justifies  the  assumption  of  a unique  effective  heat  transfer 
coefficient  for  both  the  base  and  the  lead. 

The  results  which  have  been  presented  to  this  point  were 
obtained  by  applying  the  mathematical  model  developed  for  the 
Tekform  package  used  in  the  experimental  work.  However,  other 
microelectronic  packages  which  are  commercially  available  ex- 
hibit different  characteristics.  The  present  study  applies  to 


When  h for  the  lead  and  the  base  are  different,  h..  , instead 

_ lead 

of  h is  used  in  Eq . (14)  and  hbase  instead  of  h is  used  in 
Eq.  (15). 
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-hr-  F 

Average 

o* 

max  \ 

0 

20 

40 

60 

Radial  Stress 

-0.617 

-0.650 

-0.681 

- 

-0.713 

%A  Percent  change 

0 

5.3 

10.3 

15.5 

Tangential  Stress 

0.445 

0.454 

0.460 

0.467 

%A  Percent  change 

0 

2.0 

3.3 

n 

TABLE  3 

Maximum  non-dimensional  transient  stresses  at  the  lead-glass 
interface  when  h is  different  for  the  lead  and  the  base 
compared  to  the  maximum  value  of  the  stresses  when 
R = 210  Btu/ft2-hr-«F.  (aE  - Rlead  - Rbase> 

Average  a is  calculated  by  averaging  the  values  of  the 

maximum  stresses  obtained  with  various  values  of  hlead  and 

h.  for  a given  Ah.  (e.g.,  for  Ah  = 40,  the  maximum  stresses 

D3S6  ^ 

obtained  using  h^ea£j  = 250  and  hkase  = 210,  h^ea(j  = 230  and 

h,  = 190,  h.  . = 210  and  h,  = 170  Btu/ft2-hr-°F  were 
base  lead  base 

averaged  to  obtain  the  average  cr*aj{)  . 
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all  packages  which  are  characterized  by  certain  material  property 
ratios  (Young's  modulus,  thermal  expansion  coefficients,  Poisson's 
ratio,  heat  transfer  properties)  similar  to  the  ratios  of  the  Tek- 
form  package  which  was  studied.  When  this  similarity  condition 
is  fulfilled  the  only  other  factor  that  affects  the  stress  level 
is  the  geometrical  characteristic  of  the  package.  Since  the  lead, 
the  glass  and  the  base  were  modeled  as  3 concentric  cylinders,  the 
ratios  of  the  cylinder  radii  (rQ/r^  and  r^/r 2 ) influence  the  stresses. 

The  non-dimensional  radial  and  tangential  stresses  at  the 
lead-glass  interface  were  computed  for  different  radius  ratios 
and  plotted  in  Fig.  19  and  20.  It  is  believed  that  these 
graphs  provide  the  means  for  evaluating  the  effect  of  package 
geometry  on  the  maximum  stresses  for  most  of  the  microelectronic 
packages  which  are  commercially  available. 

Fig.  19  and  20  illustrate  the  influence  of  the  choice  made 
for  the  radius  of  the  third  cylinder  which  represents  the  base 
of  the  package.  The  Tekform  package  which  was  used  in  this 
study  was  characterized  by  a ratio  t1/r2  equal  to  2.417.  From 
Fig.  19  it  can  be  determined  that  the  percent  change  in  the 
value  of  the  maximum  radial  stress  for  r^/r2  = 2.5  when  r0/r± 
is  varied  from  1.5  to  4 is  12.6  percent.  Fig.  20  shows  that 
the  percent  change  in  the  tangential  stress  corresponding  to 
the  same  value  of  r^/r2  and  for  the  same  variation  in  rQ/r^  is 
4.4  percent.  Thus,  it  appears  that  for  the  specific  package 
used  for  experimental  purposes,  the  choice  of  the  outside  radius 
of  the  third  cylinder  is  relatively  unimportant  insofar  as  the 
determination  of  the  maximum  stress  is  concerned.  It  should  be 
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Fig.  19  Maximum  non-dimensional  transient  radial  stress 
in  the  glass  at  the  lead-glass  interface  for  the 
radius  and  material  property  ratios  shown  on  the 
figure . 
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noted,  however,  that  for  other  values  of  the  corresponding  ratio 
r1/r2  it  is  possible  that  the  choice  of  rQ/r1  could  significantly 
influence  the  calculated  maximum  stress. 

The  thermal  shock  test  condition  A (32°F  to  212 °F)  of  the 
MIL- STD-88 3B  was  applied  to  the  Tekform  package  used  during  the 
experiments.  The  maximum  radial  and  tangential  stresses  occur- 
ring during  this  thermal  shock  test  were  presented  in  non-dimen- 
sional form.  The  thermal  shock  test  procedure  includes  test 
levels  which  are  more  severe  because  they  require  larger  tempera- 
ture extremes.  The  maximum  stress  values  depend  linearly  on  the 
temperature  of  the  cold  and  hot  baths.  Therefore,  the  results 
presented  for  test  A can  be  easily  calculated  for  any  other 
test  level  provided  that  the  effective  heat  transfer  coefficient 
is  the  same  as  that  determined  for  test  level  A.  The  maximum 
stress  values  corresponding  to  the  various  test  levels  for  the 
Tekform  experimental  package,  based  on  this  proviso  are  presented 
in  Table  4.  It  should  be  noted  that  5000  psi  is  a reasonable 
estimate  of  the  radial  tensile  strength  of  the  glass-lead  inter- 
face. Thus,  Table  4 indicates  that  the  higher  test  levels 

» 

could  be  damaging  to  this  interface. 
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CHAPTER  VI 

PRACTICAL  APPLICATIONS  OF  THE  RESULTS 

The  results  obtained  from  this  investigation  can  be  used 
by  the  designer  or  the  screener  of  a microelectronic  package 
to  predict  the  maximum  stresses  in  the  glass  seal  when  the  pack- 
age is  subjected  to  a thermal  shock  test. 

The  package  considered  must  have  circular  leads  and,  the 
leads  and  the  base  must  be  made  of  the  same  material.  In  order 
to  be  able  to  use  the  graphs  provided  in  the  preceding  Chapter  5 
the  package  under  study  must  have  the  same  material  property 
ratios  as  those  used  in  preparing  the  graphs  (19,20,21,22,23,24) 


IEm/Ese  ' 2-0i:  “n/“se  * 1'27:  V"se  * l-2' 
RL(Ab/N)/kse(L-L0)2  - 96.459,  P^/p^c^  - 1.681) 

The  general  procedure  which  is  to  be  followed  in  order  to 

evaluate  the  maximum  thermal  stresses  during  the  thermal  shock 


test  is  as  follows: 


1)  The  maximum  values  of  the  stresses  at  the  lead-glass 
interface  are  plotted  as  a function  of  the  radius  ratios  rQ/r^ 
and  r^/r2 • First  evaluate  r^/r2,  the  ratio  of  the  outside  radius 
of  the  glass  annulus  to  the  radius  of  the  lead.  Theoretically, 
the  outside  radius  of  the  third  cylinder  is  that  radius  for  which 


the  radial  stress  is  zero.  Determine  the  two  possible  choices 
for  r as  shown  on  Fig.  21.  Assume  [r  ].,  to  be  the  outside 
radius  of  the  annulus  representing  the  base  and  evaluate  rQ/r^ . 

2)  When  a microelectronic  package  is  manufactured,  resi- 
dual stresses  are  present  in  the  glass  seal  due  to  the  cooling 
of  the  seal  from  the  setting  point  of  the  glass  to  room  tempera- 
ture, and  the  fact  that  the  thermal  expansion  coefficients  of  the 
metal  and  the  glass  are  usually  not  equal.  Evaluate  the  non- 
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dimensional  radial  and  tangential  residual  stresses  in  the  seal 
using  Fig.  22  and  23.  Note  that  am/ase  for  this  case  is  very 
close  to  1 since  for  this  temperature  interval  the  coefficients 
of  thermal  expansion  for  the  metal  and  glass  are  matched.  Deter- 
mine the  absolute  value  of  the  stresses  using 


E o 0 

g g o 


(22) 


or  a = a *E  a 0 

g g o 


(23) 


where  0Q  represents  the  temperature  difference  THIGHER  TEMp>  BATH 
TLOWER  TEMP.  BATH* 

3)  At  the  end  of  each  thermal  shock  cycle,  after  thermal 
equilibrium  has  been  reached,  the  thermal  stresses  reach  their 
steady-state  values.  Determine  the  steady-state  non-dimensional 
radial  and  tangential  stress  when  the  package  is  taken  from  the 
lower  temperature  bath  to  the  higher  temperature  bath  by  means 
of  Fig.  24  and  25.  Compute  the  absolute  value  of  the  stresses 
using  Eq.  23. 

4)  The  first  step  in  the  thermal  shock  test  procedure  is 
to  immerse  the  package  in  the  lower  temperature  bath.  Evaluate 
the  residual  stresses  when  ten.perature  equilibrium  is  reached 

at  the  lower  temperature  bath  by  multiplying  the  absolute  values 
of  the  stresses  obtained  in  step  3 by  er/e0  where  = 

TLOW  TEMP.  BATH  " TROOM* 

5)  When  the  package  is  transferred  from  the  lower  tempera- 
ture bath  to  the  higher  temperature  bath  a transient  maximum 


stress  occurs.  Evaluate  the  transient  maximum  non-dimensional 


0N-DIMENSI0NAL  RADIAL  STRESS 

12  -0.0018  -0.0024  -0.0030  -0.0036  -0.0042  -0.0048 


Fig.  23  Non-dimensional  residual  tangential  stress  in  the 
glass  at  the  lead-glass  interface  for  the  radius 
and  material  property  ratios  shown  on  the  figure. 
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Non-dimensional  steady-state  stress  in  the  glass  at 
the  lead-glass  interface  for  the  radius  and  material 
property  ratios  shown  on  the  figure. 
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Fig.  25  Non-dimensional  steady-state  tangential  stress  in 
glass  at  the  lead-glass  interface  for  the  radius 
material  property  ratios  shown  on  the  figure. 
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radial  and  tangential  stresses  using  Fig.  18  and  19.  Compute 
the  absolute  values  of  these  stresses  with  the  help  of  Eq.  23. 

6)  The  superposition  of  the  stresses  evaluated  in  the  pre- 
ceding steps  yields  the  value  of  the  actual  radial  and  tangen- 
tial stresses  at  each  step  of  the  thermal  shock  test  procedure. 

a)  The  initial  residual  stresses  at  room  temperature  are 
equal  to  the  stresses  obtained  in  step  (2) . 

b)  The  residual  stresses  at  the  lower  temperature  bath 
are  equal  to  the  stresses  obtained  in  step  (2)  + the  stresses 
obtained  in  step  (4) . 

c)  The  maximum  transient  stresses  occurring  when  the  pack- 
age is  transferred  from  the  lower  temperature  bath  to  the  higher 
temperature  bath  are  equal  to  the  stresses  obtained  in  step  (6-b) 
+ the  stresses  obtained  in  step  (5). 

d)  The  steady-state  stresses  at  the  higher  temperature  bath 
are  equal  to  the  stresses  obtained  in  step  (3)  + the  stresses 
obtained  from  step  (6-b) . 

e)  The  maximum  transient  stresses  occurring  when  the  pack- 
age is  transferred  back  in  the  lower  temperature  bath  are  equal 
to  the  stresses  obtained  from  step  (6-d)  + (-the  stresses  ob- 
tained from  step  (5)^). 

7)  It  was  stated  in  step  1 that  the  third  outside  radius 


When  the  absolute  value  of  the  stresses  is  obtained  in  step  5 

eo  ls  THIGH  TEMP.  BATH  “ TLOW  TEMP.  BATH*  In  Step  (6_c)  0O 
becomes  TLQW  TEMp>  BATH  - THIGH  TEMPi  BATH  which  is  the  negative 


value  of  eQ  used  in  step  5. 


Therefore  the  absolute  value  of 
the  maximum  transient  stresses  when  the  package  is  transferred 
from  the  higher  to  the  lower  temperature  bath  is  the  negative 
of  the  ones  obtained  in  step  5. 
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representing  the  base  can  have  two  values.  To  compare  the 
effects  of  the  two  possible  radii  on  the  thermal  stress  history 
of  the  seal,  replace  [rQ]^  by  [rQ]  and  repeat  the  entire  pro- 


cedure until  the  end  of  step  6. 

Repeating  the  transfer  of  the  package  from  the  lower  temper- 
ature bath  to  the  higher  temperature  bath  and  back  again  as  re- 
quired by  the  thermal  shock  test  procedure  yields  the  maximum 
transient  stresses  determined  in  steps  (6-c)  and  (6-e) . 

The  following  example  illustrates  the  procedure  explained 
above . 

EXAMPLE:  An  Isotronics  flatpack  type  package  with  the  lead- 

glass-base  configuration  as  shown  in  Fig.  26  is  subjected  to 
the  test  condition  B of  the  thermal  shock  test.  The  temperature 
extremes  for  this  test  are  given  as  -67°F  and  257°F.  The  leads 
and  the  base  are  made  of  Kovar  and  it  is  assumed  that  7052  Corning 
glass  is  used  for  sealing  the  leads  to  the  base.  Using  the 
procedure  described  above,  the  history  of  the  maximum  thermal 
stresses  is  calculated. 

Step  1 - The  radius  of  the  lead  is  0.01"  and  the  radius  of  the 
glass  is  0.035".  Therefore  the  ratio  is 


r1/r2  ==  3.5 


The  reasonable  candidates  for  the  assumed  outside  radii  are 


[r  ] = .06"  and  [r  ] = .05".  [r  ] is  chosen  first  and  the 

L oJ^  L oJ2  o i 


other  radius  ratio  becomes 

UoVh  ' 1,71 


1 


i 


Fig.  26  Schematic  Representation 
through  and  the  base  of  a 
Package . 
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Step  2 - The  non-dimensional  residual  stresses  after  cooling  to 

room  temperature  are 

(Using  Fig.  22)  Radial  stress  = -.00146 
(Using  Fig.  23)  Tangential  stress  = -0.0139 

The  materials  (Kovar  and  glass)  being  the  same  as  the  ones 

used  in  this  study,  the  ratios  of  material  properties  are 

equal  to  those  on  Fig.  22  and  23.  The  absolute  value  of  the 

stresses  is  evaluated  using  Eq.  (23)  where  Eg  = 8.2  x 106  psi, 

a = 2.95  x 10-6/°F,  0 = 77  - 817  = -740°F  for  this  step, 

g o 

The  absolute  values  for  the  stresses  are 

Radial  stress  = -.00146  x 8.2  x 10^  x 2.95  x 10  ^ x 
-740  - 26.13  psi 

Tangential  stress  = -0.0139  x 8.2  x 10^  x 2.95  x 10  ^ x 

-740  = 248 .82  psi. 

Step  3 - The  non-dimensional  steady-state  stresses  at  the  higher 
temperature  bath  are 

_2 

(Using  Fig.  24)  Radial  stress  = 3.87  x 10 

-2 

(Using  Fib.  25)  Tangential  stress  = 36.3  x 10 

The  absolute  values  of  these  stresses  using  E^  = 8.2  x 10^  psi, 

a = 2.56  x 10_6/°F  and  6 = 257  - (-67)  = 324°F  is 

g o 

Radial  stress  = 3.87  x 10  ^ x 8.2  x 10^  x 2.56  x 10  ^ x 

324  = 263.21  psi 

Tangential  stress  = 36.3  x 10~  x 8.2  x 10  x 2.56  x 10  x 

324  = 2468.91  psi. 

Step  4 - The  steady-state  stresses  after  immersing  the  package 

from  room  temperature  in  the  low  temperature  bath  are  evaluated 
using 

0 = -67  - 77  = 144°F 

r 
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Thus  using  the  absolute  values  of  the  stresses  obtained  in 


Step  3 these  stresses  are 


-144 

Radial  stress  = 263.21  x . = -116.98  psi 


Tangential  stress  = 2468.91  x — = -1097.29  psi. 


Step  5 - The  non-dimensional  maximum  transient  thermal  stresses 


(Using  Fig.  19)  Radial  stress  = -55.7  x 10 


(Using  Fig.  20)  Tangential  stress  = 52.4  x 10 


with  absolute  values 


Radial  stress  = -55.7  x 10  ^ x 8.2  x 10^  x 2.56  x 10  ^ x 


324  = -3788.38  psi 


— 2 6 —ft 

Tangential  stress  = 52.4  x 10  x 8.2  x 10  x 2.56  x 10  x 


324  = 3563.94  psi. 


Step  6 - The  sequence  of  the  thermal  stress  is : 


a)  Initial  residual  stresses: 


= 26.13  psi 


a = 248.82  psi 


b)  Steady-state  stresses  at  low  temperature  bath: 


ar  = 26.13  - 116.98  = -90.85  psi 


o = 248.82  - 1097.29  = -848.47  psi 


c)  Maximum  transient  stresses  going  from  low  to  high 


temperature  bath 


o r = -90.85-3788.38  = -3879.23  psi 


(7  = -848.47  + 3563.94  = 2715.47  psi 


d)  Steady-state  stresses  at  higher  temperature  bath: 


j = 263.21  - 90.85  = 172.36  psi 


o = 2468.91  - 848.47  = 1620.44  psi 


=L 
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e)  Maximum  transient  stresses  going  from  high  to 
low  temperature  bath: 

or  = 172.36  - (-3788.38)  = 3960.74  psi 

a = 1620.44  - 3563.94  = -1943.50  psi 

<t> 

This  cycle  is  shown  in  Fig.  26. 

Step  7 - The  same  procedure  is  repeated  using  [r  ] . In  this 

o 2 

case  [r  ] /r-,  = 1.43.  Since  this  ratio  is  not  given  in  the 
° 2 1 

graphs  it  will  be  assumed  that  rQ/r^  = 1.5.  ri/r2  s 
same  as  before,  i.e.,  3.5. 

Step  2 - The  non-dimensional  residual  stresses  at  room  tempera- 
ture are: 

(Using  Fig.  22)  Radial  stress  = -.00040 

(Using  Fig.  23)  Tangential  stress  = -.134  x 10  1 

The  absolute  values  are 

Radial  stress  = -.00040  x 8.2  x 10^  x 2.95  x 10  ® x 
-740  = 7.16  psi 

Tangential  stress  = -.134  x 10  ^ x 8.2  x 10^  x 2.95  x 10  ® x 

-740  = 239.87  psi 

Step  3 - The  non-dimensional  steady-state  stresses  at  the  higher 

temperature  bath  are: 

(Using  Fig.  24)  Radial  stress  = .0107 
(Using  Fig.  25)  Tangential  stress  = .355 

The  absolute  values  of  these  stresses  are 

Radial  stress  = .0107  x 8.2  x 10^  x 2.56  x 10  ® x 

324  = 72.78  psi 

Tangential  stresses  = .355  x 8.2  x 106  x 2.56  x 10  ^ x 

324  = 2414.50  psi 

Step  4 - The  steady-state  stresses  after  immersing  the  package 
from  room  temperature  in  the  low  temperature  bath  are 
evaluated  using 


Radial  Stress 


Thermal  stresses  for  test  condition  B (-67 °F 
when  r = [r  1 = 0.06". 


A = stress  at  room  temperature;  B = stress  at  lower 
temperature  bath;  C = maximum  transient  stress  when 
the  package  is  immersed  in  the  higher  temperature 
bath;  D = maximum  transient  stress  when  the  package 
is  immersed  back  in  the  lower  temperature  bath; 

C'  = repetition  of  low-to-high  cycle;  D'  = repetition 
of  high-to-low  cycle. 


i 


! 

I 


1 
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6 = -67  - 77  = -144°F 

r 

Thus  using  the  absolute  values  of  the  stresses  obtained  in 
Step  3 these  stresses  are 


Radial  stress  = 72.78  x = - 32.35  psi 

-144 

Tangential  stress  = 2414.50  x = “1073. 11  psi 

Step  5 - The  non-dimensional  maximum  transient  thermal  stresses 
are: 

(Using  Fig.  19)  Radial  stress  = -.547 


(Using  Fig.  20)  Tangential  stress  = .526 


The  absolute  values  of  these  stresses  are 


Radial  stress  = -.547  x 8.2  x 10^  x 2.56  x 10  ^ x 
324  = -3720.37  psi 

C.  __  r 

Tangential  stress  = .526  x 8.2  x 10  x 2.56  x 10  x 

324  = 3577.54  psi 

Step  6 - The  sequence  of  the  thermal  stresses  is: 

a)  Initial  residual  stresses 

0^  = 7.16  psi 

a,  = 239.87  psi 

(j> 

b)  Steady-state  stresses  at  low  temperature  bath: 

ar  = 7.16  - 32.35  = -25.19  psi 

a , = 239.87  - 1073.11  = -833.24  psi 

c)  Maximum  transient  stresses  going  from  low  to  high 
temperature  bath: 

o r = -25.19  - 3720.37  = -3745.56  psi 
oA  = -833.24  + 3577.54  = 2744.30  psi 

d)  Steady-state  stresses  at  higher  temperature  bath: 

or  = 72.78  - 25.19  = 47.59  psi 
0 = 2414.50  - 833.24  = 1581.26  psi 

e)  Maximum  transient  stresses  going  from  high  to  low 
temperature  bath: 

<t>r  = 47.59  - (-3720.37)  = 3767.96  psi 
o.  = 1581.26  - 3577.54  = -1996.28  psi 


This  cycle  is  shown  in  Fig.  27. 


Thermal  stresses  for  test  condition  B (-67°F 
when  r =[r]  =0.05". 


A = stress  at  room  temperature;  B = stress  at  lower 
temperature  bath;  C = maximum  transient  stress  when 
the  package  is  immersed  in  the  higher  temperature 
bath;  D = maximum  transient  stress  when  the  package 
is  immersed  bark  in  the  lower  temperature  bath; 

C'  = repetition  of  low-to-high  cycle;  D'  = repetition 
of  high-to-lov:  cycle. 
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CHAPTER  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 
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The  mathematical  model  developed  in  this  study  enables  one 
to  evaluate  the  transient  thermal  stresses  that  are  experienced 
in  the  annular  glass  seal  of  a microelectronic  package  with 
circular  leads  when  it  is  subjected  to  a thermal  shock  test. 

The  results  show  that  the  transient  radial  stresses  occurring 
during  the  thermal  shock  test  which  are  caused  by  large  tempera- 
ture gradients  can  become  significantly  higher  than  the  stresses 
associated  with  a uniform  temperature  distribution  which  are 
due  to  a mis-match  in  the  thermal  expansion  coefficients  of  the 
metal  and  the  glass.  When  the  transient  tangential  stress  is 
compared  with  the  stresses  that  would  occur  due  to  only  a mis- 
match in  thermal  expansion  coefficients  and  &se,  the  difference 
is  shown  to  be  very  small.  The  comparison  made  above  can  be 
illustrated  using  the  results  of  the  example  given  in  the  pre- 
ceding chapter.  The  transient  maximum  radial  stress  when  the 
package  is  heating  was  calculated  to  be  -3788.38  psi  whereas 
the  steady-state  stress  value  at  thermal  equilibrium  was 
263.21  psi.  For  the  tangential  stresses,  the  maximum  tran- 
sient stress  was  3393.90  psi  while  the  steady-state  value  was 
2468.91  psi. 

Information  concerning  the  expected  values  of  the  maximum 
transient  thermal  stresses  as  influenced  by  the  geometry  of  the 
seal,  the  effective  heat  transfer  coefficient  and  the  different 
thermal  shock  test  conditions  as  presented  in  this  study  is  im- 
portant in  two  ways.  From  a design  point  of  view,  the  factors 
affecting  the  stress  level  can  be  adjusted  in  order  to  reduce 
the  maximum  stresses  so  that  the  possibility  of  failure  of  the 
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seal  could  be  diminished.  From  the  screener's  point  of  view, 
the  results  can  be  used  to  predict  what  regions  of  the  seal 
will  be  stressed  significantly  and  how  severe  the  test  will  be 
when  a given  thermal  shock  procedure  is  applied  to  the  package. 
Conversely,  when  a certain  stress  level  is  desired  in  the  seal, 
the  proper  test  level  can  be  selected  by  the  use  of  the  results 
of  the  present  study. 

The  APL  programs  which  have  been  developed  to  compute 
the  transient  thermal  stresses  have  been  described  in  Chapter  V. 

The  inputs  for  these  programs  (consisting  of  the  Young's  moduli 
for  the  3 materials,  the  coefficients  of  thermal  expansion  for 
the  3 materials,  the  Poisson's  ratios  for  the  3 materials,  the 
3 radii  corresponding  to  the  3 cylinders,  the  densities  for  the 
3 materials,  the  specific  heats  for  the  3 materials,  the  thermal 
conductivity  of  the  sealing  material,  the  outside  surface  area 
of  the  lead  in  contact  with  the  fluid,  the  outside  surface  area 
of  the  base  and  lid,  the  volume  of  one  lead,  the  total  volume  of 
the  base  and  lid,  the  initial  temperature  of  the  package,  the 
final  temperature  of  the  package,  the  effective  heat  transfer 
coefficient,  the  radius  where  the  stresses  are  to  be  evaluated, 
the  total  number  of  leads,  the  time  increment  and  the  total  num- 
ber of  time  increments)  are  not  restricted  to  materials  with  pro- 
perty ratios  equal  to  those  specified  on  the  graphs  presented  in 
the  preceding  chapters.  The  computer  programs  may  be  used  to 
determine  the  thermal  stress  history  of  any  package  during  a 
thermal  shock  test  provided  that  the  leads  are  circular.  An 
example  illustrating  the  use  of  the  program  is  given  in  Appendix  C. 

The  work  herein  reported  provides  a method  for  estimating 
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the  stresses  which  are  present  in  the  vicinity  of  a lead-through 
seal  of  a microelectronic  package  when  subjected  to  a thermal 
shock.  The  following  recommendations  are  made  for  further  work 
with  the  objectives  of  verifying  the  validity  of  certain  sim- 
plifying assumptions  which  were  made  in  the  development  of  the 
mathematical  model  and  extending  the  work  to  include  other  im- 
portant aspects  of  the  general  problem  of  estimating  the  con- 
ditions which  would  cause  failure  of  the  lead-through  seal  during 
thermal  shock. 

It  is  suspected  that  failure  of  the  glass-metal  seal  will 
occur  when  the  tensile  radial  stress  at  the  glass-metal  inter- 
face reaches  a critical  value.  At  the  present  time  little  is 
known  about  the  magnitude  of  this  critical  value  and  how  it  is 
influenced  by  the  manufacturing  process.  It  is  recommended  that 
experimental  studies  be  carried  out  to  determine  the  strength 
of  the  seal.  Subsequently,  experiments  should  be  performed  to 
verify  the  predictions  of  the  present  model  with  regard  to  seal 
failure.  This  could  be  accomplished  by  testing  a variety  of 
different  microelectronic  packages  at  different  thermal  shock 
levels . 

In  order  to  investigate  the  validity  of  the  assumptions  in- 
herent in  the  present  model  it  is  suggested  that  a more  refined 
finite  element  analysis  be  performed  for  a typical  situation 
covered  by  the  present  theory.  It  would  be  especially  interesting 
to  know  whether  the  assumptions  which  have  been  made  with  regard 
to  neglecting  the  axial  temperature  distribution  in  the  glass 
lead  to  serious  error.  A further  important  assumption  made  in 
the  mathematical  model  was  the  modeling  of  the  base  as  a third 
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outside  cylinder,  having  a certain  assumed  outer  radius.  A 
more  realistic  evaluation  of  the  stresses  could  be  obtained  if 
a finite  element  analysis  were  made  to  evaluate  the  thermal 
stresses  in  the  seal  taking  into  account  the  real  shape  of  the 

package  including  the  influence  of  any  neighboring  seals.  This  j 

would  provide  an  evaluation  of  the  accuracy  of  the  3-cylinder 
model  used  in  the  present  study. 

On  the  experimental  side,  the  uncertainties  involved  in  the 
determination  of  the  effective  heat  transfer  coefficient  for  the 
package  have  been  discussed  in  Chapter  IV.  A more  refined  ex- 
perimental technique  incorporating  the  effects  of  other  package 
variables  such  as  the  package  contents  should  be  undertaken  in 
order  to  provide  a more  accurate  determination  of  the  true  effec- 
tive heat  transfer  coefficient  for  test  condition  A of  the  MIL- 
STD-883B.  Improved  results- may  be  achieved  by  assuming  that 
the  effective  heat  transfer  coefficient  is  a function  of  the 
temperature  difference  between  the  fluid  and  the  package. 

Investigation  should  be  carried  out  for  other  test  levels 
described  in  MIL-STD-883B  to  determine  the  effective  heat  trans- 

3 

fer  coefficient  for  the  other  fluids  during  the  specified  test 
conditions . 
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In  view  of  the  fact  that  many  microelectronic  packages  have 
leads  which  are  not  circular,  further  analysis  should  be  performed 
to  provide  means  for  estimating  thermal  stresses  for  the  case 
of  rectangular  leads  or  other  lead-through  geometries. 

It  was  shown  in  the  practical  applications  of  the  results 
that  the  magnitude  of  the  thermal  stresses  in  the  seal  depended 
on  the  initial  residual  stresses  in  the  seal  caused  by  the 
cooling  of  the  package  to  room  temperature  after  the  sealing 
process  is  performed  at  high  temperatures.  In  view  of  the 
fact  that  these  residual  stresses  affect  directly  the  maximum 
thermal  stresses,  further  investigation  is  needed  to  evaluate 
accurately  the  initial  residual  stresses  in  the  glass-to- 
metal  seal. 


APPENDIX  A 


FINITE  ELEMENT  ANALYSIS  FOR  EVALUATING 
THE  TRANSIENT  TEMPERATURE  DISTRIBUTION 
IN  THE  GLASS  AND  RELATED  APL  PROGRAMS 


A.l  Evaluation  of  the  Transient  Temperature  Distribution  in 
the  Glass  Using  the  Finite  Element  Analysis 

A finite  element  analysis  is  used  to  solve  the  partial  differ- 
ential equation  (12)  in  Chapter  III  together  with  the  boundary 
conditions  (14)  and  (15).  The  mathematical  development  of  the 
analysis  can  be  found  in  reference  (11). 

The  finite  element  model  used  consisted  of  five  circular 
elements  dividing  the  glass  seal.  (See  Fig.  28) . This  model 
requires  that  6 nodal  temperatures  be  evaluated  using  the  finite 
element  method.  T^  and  Tg  correspond  respectively  to  the  tempera- 
tures of  the  lead  and  the  base.  The  intermediate  temperatures 
determine  the  radial  temperature  distribution  through  the  glass 
seal.  This  set  of  temperatures  can  be  represented  by  a matrix. 


The  mathematical  development  of  the  analysis  applied  to  the 
partial  differential  equation  (12)  leads  then  to  the  following 
system  of  equations 

[C]  + [K]{T>  + [f]  = 0 (Al) 

d T 

where  [c]  is  the  total  "capacitance"  matrix  expressed  for  each 
element  as 


1 


Fig.  28.  Model  Used  for  Finite  Element  Analysis 
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[c]  = 


. _ (2r_5-20r_2r1 3+30ror. 

: n >'■  C 2 2 1 2 1 

[r,-r. ) * . . . 

2 1 |_  symmetrical 


1“-12r15) (3r25-5r2 


(12r25-30r2‘4 


“rl+5r2rl4'3rl5] 

r1+20r23r12-2r15J 


and  r2  being  respectively  the  inner  and  outer  radius  of  each 

element.  For  the  inner  boundary  element,  [c]  = p^Cj^V^  Ijjj  is 

added  to  [c]  above.  For  the  outer  boundary  element,  [c]  = 
p,c,,V,  n n 

K - [„  :]  is  added  to  [c]  above.  [ K]  is  the  total  "conductivity1 

N 0 1 

matrix  and  for  each  element 


_ 2*k(r23-r13)  r 1 -ri 

LkJ  - -llr2-r1)^-  |_-1  l] 


For  the  inner  boundary  element,  [k]  = hAQ  [J  is  added  to  fk] 
above.  For  the  outer  boundary  element,  [k]  = [°  °]  is  added 

to  [k]  above.  [F]  is  the  "forcing"  matrix  and  applies  to  the 
boundary  elements.  For  the  inner  boundary 


[f]  = -hAoTj3] 

hAb  0 

For  the  outer  boundary  [f]  = - -jp  T^t^] 

Equation  (Al)  is  solved  in  the  time  domain  using  a finite 
difference  method. 

The  time  derivative  is  represented  by: 

Hrr  = it  <{T}i  - {T}o>  {A2) 

where  At  is  a finite  time  increment  and  f T } ^ and  {T}Q  are  the 
model  temperature  matrices  corresponding  to  the  beginning  and  end 
of  the  time  interval  At.  At  the  mid  point  of  the  time  interval 
the  temperature  matrix  can  be  written  as 


(A3) 


(T)  = ±({T}.  + (T)  ) 

Z 1 o 

Replacing  (A2)  and  (A3)  in  (Al)  and  rearranging  terms,  one  obtains 
<[k]  + ^ [C]){T)1  = (|7[C]-[k]){T}o  - 2[F]  (A4) 

Equation  (A4)  is  solved  to  yield  the  nodal  temperature  values 
at  time  t + At  (represented  by  {T}^)  given  the  nodal  values 
at  time  t (represented  by  (T}q) . The  initial  temperature 
distribution  is  used  here  to  start  the  evaluation  process. 

The  calculations  based  on  this  analysis  were  carried  out 
using  the  APL  programs  listed  in  Sections  A. 2 and  A. 3. 


A. 2 Programs  "TEMPDIS"  and  "FEM"  to  evaluate  the  [K],[C]  and 
[F]  matrices. 


VTCMP-DXSCOlV 

W w temfbxs  dt;m;mc;cm;i.oc;fow}Cou;  k ;c;Fc;TXN;ME>;M;n;ni  A3>A4>r>rp'l  ?F'2r  tf ;ta 
C13  *«-6 

1 23  »'E*-5 

133  RR«-RC33*P'C33-*’<<*:C23-*C33)  +Me>* < l*E> 

C43  HOf-(He,2)rO 
C53  **1 

C63  MNH+I 
C 73  MOPCin*-*»J 

C83  -*'-lxi<x».x  + i)iNE 

C93  *r«-<N,l)fO 

C 103  CTfKTf(ll,H)rO 
Cl  1 3 **-«-l 

C123  eCi-1,5 

Cl 33  i>2:»"cm  »t 
C 143  »-5:«-oc«.2r0 

C15J  coccvSl^-f'ooC^ll 

C 163  H-0 

C 173  *-31  x«-*  + l 
C183  ROWt.LOC[i] 

C 193  J«-0 

C203 

C213  eoi.4-LoecJ3 

C223  KT[F:OW}COl.3«.ICT[ROIi»  JCOU3  + «C  * > J3 
C233  CT[Fow;eoi.3«.CTCF:ow;coi.3+ccx  > J3 
C243  -*l-4x  t < J<2) 

C253  -*«.7x»el«ec 

C263  •♦•-8 

1273  *-7SFCROW»  1 3*-rC*owH3+F*CI > 13 
£283  »-S;-»‘-3*l(*«2) 

C 293  ec^ee+ j 
C303  -»i-2x|Ei.<nc 

C313  "♦-•‘T 

C323  c«-ct 
w 


2 r l 


VFCMCQ3V 

V FEM  DTfXjKV|«f jCViCt 

C13  R«-RR£MOP£El.J33 

•C23  «v«.(<o2)xXftx<<FC23*3)-(«C13*3))+3x(»*C23-*C13>»2)x(2 

C33  «*«-  2 2 tO 

143  EfKVf KS 
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A. 3 Program  "GTS"  to  compute  the  transient  temperature  distribu' 
tion  in  the  glass  seal. 
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B.l  Main  program  "STR"  to  compute  the  transient  thermal  stresses 
at  a given  radius. 
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B.2  Programs  "STRESS 


and  "Z 


to  compute  the  thermal  stresses 


at  each  increment  of  time. 
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APPENDIX  C 

APL  FUNCTION  EXPLAINING  HOW  THERMAL  STRESSES  ARE 
COMPUTED  USING  THE  APL  PROGRAMS  PRESENTED  IN 
APPENDICES  A AND  B 
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